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Technical  Section 


Technical  Objectives 

The  objective  of  the  work  was  to  investigate  nanoengineer ed  surfaces  that  could  modulate  the  wetting  and 
phase  change  processes  to  enhance  heat  transfer  performance.  Phase  change  processes,  associated  with 
the  large  latent  heat,  have  been  widely  used  in  applications  where  high  heat  fluxes  and  high  heat  transfer 
coefficients  are  desired.  One  example  is  thermal  management  which  has  become  a  critical  bottleneck  for 
the  advancement  of  a  variety  of  important  defense,  space,  and  commercial  systems.  Current  state-of-the- 
art  cooling  approaches  have  significant  limitations,  particularly  for  high  heat  flux  applications.  Two- 
phase  fluidic  cooling  systems,  utilizing  liquid-vapor  phase-change  absorbs  large  fluxes  with  minimal 
changes  in  device  temperature  [1-3],  has  been  considered  as  one  of  the  most  promising  approach  to 
address  the  demands  for  thermal  management.  While  efforts  have  focused  primarily  on  boiling  flows  in 
microchannels,  liquid- vapor  instabilities  during  phase-change  in  these  systems  lead  to  local  dry-out,  non- 
uniform  temperature  distributions,  and  significant  decreases  in  critical  heat  flux  [4-8].  Despite  extensive 
research,  there  remain  significant  hurdles  in  the  implementation  of  this  technology.  In  contrast,  two- 
phase  jet  impingement  and  spray  cooling  techniques  have  been  investigated  to  a  much  lesser  extent,  even 
though  such  a  methods  promise  high  heat  transfer  coefficients  (>100  W/cm2K)  with  high  heat  dissipation 
capability  (>1000  W  and  >1000  W/cm2)  via  thin  film  evaporation  [1,  9].  The  implementation  has 
comparable  challenges.  Previous  experimental  attempts  have  typically  led  to  pool  boiling  due  to  chamber 
flooding,  or  liquid  dry-out  due  to  insufficient  liquid  supply  [1,  10].  These  undesired  effects  significantly 
decreased  heat  transfer  coefficients  and  heat  removal  rates.  In  one  aspect  of  this  work,  we  took  advantage 
of  state-of-the-art  nanoengineering  capabilities  to  address  these  challenges.  We  obtained  fundamental 
understanding  of  the  complex  fluidic  and  thermal  transport  processes  on  nanoengineered  surfaces  to 
realize  the  great  potential  of  thin  film  evaporation. 

Other  than  liquid- vapor  phase-change,  we  also  explored  modification  of  surface  wettability  to  enhance 
condensation  heat  transfer.  Dropwise  condensation,  where  the  condensate  forms  discrete  droplets  rather 
than  continuous  films  covering  the  substrate,  is  considered  as  one  of  the  most  effective  approaches  to 
enhance  the  heat  transfer  coefficient.  Previous  work  has  demonstrated  the  application  of  nanostructured 
superhydrophobic  surfaces  where  condensate  can  be  spontaneously  removed  via  a  surface-tension- driven 
mechanism  [11].  However,  the  nucleation  density  on  these  surfaces  is  relatively  low  since  the  phase 
change  process  relies  on  high  energy  active  sites  to  initiate  nucleation  at  low  super  saturations,  limiting  the 
overall  heat  transfer  performance.  Furthermore,  air  pockets  trapped  beneath  the  droplets  during  growth 
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reduce  the  contact  area  between  the  condensing  droplet  and  substrate,  which  increases  the  thermal 
resistance  and  reduces  the  heat  transfer  coefficient  [12].  In  this  work,  we  demonstrated  that  by  infusing 
hydrophobic  oil  into  superhydrophobic  surfaces,  the  nucleation  density  was  significantly  increased  while 
maintaining  easy  droplet  removal.  Such  a  surface  could  potentially  be  the  ideal  surface  for  condensation 
heat  transfer. 

Technical  Approach 

Fig.  1  shows  the  envisioned  hybrid  nanostructured  surface  to  accommodate  high  heat  flux  applications 
with  local  flux  nonuniformities.  Based  on  our  preliminary  findings  [13,  14],  the  features  include: 

•  Pillared  structures  (Fig.  la)  that  serve  as  capillary  wicks  where  the  height,  spacing,  and 
distribution  control  the  rate  of  liquid  propagation; 

•  Asymmetric  structures  (Fig.  lb)  that  can  maneuver  the  direction  of  liquid  propagation.  In  the 
presence  of  a  hotspot  where  the  local  flux  is  higher,  the  structures  direct  more  liquid  towards  the 
region  for  enhanced  heat  removal; 

•  Structures  with  fine  features  (Fig.  lc)  to  allow  three-dimensional  control  and  variations  in  liquid 
film  thickness. 


jet 


Fig.  1  Schematic  of  envisioned  hybrid  nanostructured  surface  to  accommodate  high  heat  flux  applications  with  local 
flux  nonuniformities  including  a)  pillared  structures,  b)  asymmetric  structures,  and  c)  structures  with  fine  features. 


A  key  step  to  the  realization  of  the  functional  surface  envisioned  in  Fig.  1,  however,  is  a  systematic 
investigation  and  fundamental  understanding  of  the  fluid-structure-heat  interaction  and  rich  interfacial 
phenomena  within  these  complex  configurations.  The  goal  of  our  work  was  to  fabricate  well-controlled 
and  new  nanostructure  geometries,  and  pursue  a  series  of  experiments  and  modeling  work  to  determine 
how  various  parameters  (shape,  size,  and  spacing,  etc.)  of  the  nanostructured  features  affect  the  coupled 
heat  and  mass  transport  processes.  The  understanding  gained  from  this  work  combined  with  the 
parametric  studies  provided  new  functionalities  and  offered  important  steps  towards  the  fine  control 
necessary  to  achieve  the  highest  cooling  capability  with  thin  film  evaporation. 

Summary  of  Achievements 

Through  the  duration  of  the  research  grant,  we  have  made  several  contributions: 

1.  We  developed  a  predictive  model  for  liquid  propagation  rates  based  on  the  diameter,  height,  and 
period  of  micropillar  arrays,  which  is  validated  with  experiments.  This  model  offers  an  important 
step  towards  optimizing  pillar  designs  to  achieve  high  flux  thin  film  evaporation. 

2.  We  demonstrated  that  the  liquid  droplets  could  exhibit  non- wetting  behavior  on  superhydrophilic 
microstructured  surfaces  at  superheated  conditions.  Vapor  is  generated  by  surface  heating  and 
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accumulates  beneath  the  droplet  before  flowing  through  the  pillar  arrays  and  the  increased 
pressure  of  vapor  supports  the  droplet  on  top  of  silicon  micropillar  arrays  in  a  non- wetting  state. 
Such  non- wetting  state  is  unfavorable  for  thin  film  evaporation.  A  force  balance  model  was 
developed  to  delineate  the  wetting  from  the  non- wetting  regime  of  droplets  based  on  the  droplet 
size  and  the  pillar-array  geometry.  This  study  is  important  in  outlining  conditions  to  avoid  non¬ 
wetting  of  liquids  in  thin  film  evaporation. 

3.  We  developed  a  new  metrology  technique,  Pulsed  Evaporative  Transient  Thermometry  (PETT), 
to  obtain  both  the  instantaneous  heat  capacity  and  heat  transfer  conductance,  i.e.,  the  product  of 
the  heat  transfer  coefficient  and  surface  area,  during  liquid  evaporation  from  micro  and 
nanostructures.  This  approach  experimentally  captures  the  transient  thermal  response  as  the  liquid 
evaporates  from  micro  and  nanostructures  by  a  series  of  applied  heat  pulses,  and  uses  a  lumped 
system  model  to  extract  these  two  quantities.  We  demonstrated  the  applicability  of  PETT  to 
capture  the  evaporation  process  for  two  model  systems,  copper  microwires  and  nanoporous 
alumina  membranes  with  a  coated  metal  film.  This  PETT  technique  promises  to  be  a  valuable 
quantitative  approach  to  fundamentally  investigate  the  relationship  between  wetting  morphology 
and  heat  transfer  behavior  for  a  variety  of  micro  and  nanostructured  devices. 

4.  We  demonstrated  using  nanoporous  membranes  to  achieve  thin  film  evaporation  over  large  areas. 
Steady  and  transient  heat  transfer  coefficient  over  104  W/m2K  was  obtained  with  isopropyl 
alcohol  (IP A)  as  the  coolant,  which  is  higher  than  that  of  pool  boiling  and  is  less  than  two  orders 
of  magnitude  smaller  than  the  theoretical  limit.  This  work  offers  insights  into  optimal 
experimental  designs  towards  achieving  kinetic  limits  of  heat  transfer  for  thin  film  evaporation 
based  thermal  management  solutions. 

5.  We  showed  that  the  nucleation  density  on  oil-infused,  silane-coated  structured  surfaces  can  be 
significantly  increased  by  the  use  of  disordered  long-chain  silane  coatings  that  result  in  nucleation 
sites  limited  only  by  the  density  of  pillar  structures  comprising  the  surface.  We  find  that  the 
increase  in  the  nucleation  density  can  be  explained  by  heterogeneity  in  the  surface  energy  of  the 
silane  coating  and  the  reduced  water-oil  interfacial  energy.  This  effect  could  potentially  be  used 
to  significantly  improve  the  heat  transfer  coefficient  in  condensation  by  controlling  the  nucleation 
density. 

Detailed  contributions 

1)  Liquid  Propagation  on  Micropillar  Arrays 

Nanostructures  with  novel  functionalities  provide  unique  opportunities  to  achieve  the  desired  fluidic 
control  to  enhance  heat  transfer  performance.  We  have  successfully  fabricated  nanostructured  surfaces 
that  can  manipulate  liquids  in  desired  directionalities.  However,  to  successfully  design  surfaces  for 
optimized  heat  transfer,  we  need  to  obtain  quantitative  understanding  of  liquid  dynamics  on  engineered 
surfaces  that  can  maneuver  the  liquid  and  provide  the  necessary  liquid  flow  rates  to  avoid  dry -out.  We 
have  initially  focused  our  studies  on  uniformly  spaced  micropillar  arrays,  which  we  determined  to  be 
complex  and  rich  of  interfacial  phenomena.  We  developed  a  predictive  model  for  liquid  propagation  rates 
based  on  the  diameter,  height,  and  period  of  micropillar  arrays,  which  is  validated  with  experiments.  This 
model  offers  an  important  step  towards  optimizing  pillar  designs  to  achieve  high  flux  thin  film 
evaporation. 
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We  have  fabricated  a  range  of  micropillar  arrays  with  various  geometries.  The  pillar  diameters  ranged 
from  d  -  1.3  -  6.4  pm,  the  periods  ranged  from  /  =  4.5  -  30  pm,  and  the  heights  ranged  from  h-  8.3  - 
20  pm.  An  example  of  a  uniform  square  array  of  cylindrical  pillars  with  diameter  d  -  2.3  pm,  period 
(pillar  center-center  distance)  /  =  4.5  pm,  and  height  h  -  8.3  pm,  is  shown  in  Fig.  2(a).  The  pillars  were 
defined  by  photolithography  and  etched  by  deep  reactive  ion  etching  (DRIE)  in  silicon.  Due  to  the 
alternating  etch  and  passivation  steps  characteristic  of  DRIE,  scallop  features  were  formed  on  the  sides  of 
the  pillars  (Fig.  2  (b)).  


(a)  (b) 


Fig.  2  Scanning  electron  micrograph  (SEM)  of  a  representative  micropillar  array,  (a)  Micropillars  with  diameter 
d  =  2.3  pm,  period  /  =  4.5  pm,  and  height  h  =  8.3  pm.  (b)  Magnified  image  of  a  single  micropillar  from  (a)  showing 
scallop  features  on  the  sides  as  a  result  of  deep  reactive  ion  etching  (DRIE). 


The  developed  model  determines  the  liquid  propagation  rate  in  the  pillar  arrays  based  on  the  balance 
between  the  capillary  pressure  and  the  viscous  resistance  following  Washburn  dynamics.  We  determined 
the  capillary  pressure  in  the  micropillar  array  using  a  thermodynamic  approach  minimizing  interfacial 
free  energy.  The  pressure  is  defined  as  the  change  in  surface  energy  per  unit  volume, 

Pcap=AE/AV,  (1) 


(a)  (b)  (c) 


Fig.  3.  Comparison  of  the  water  meniscus  shape  from  simulation  and  experiment,  (a)  Simulated  meniscus  shape 
using  Surface  Evolver  (SE)  where  d  =  6.4  pm  and  /  =  15  pm.  (b)  Image  of  the  water  meniscus  obtained  using 
interference  microscopy  in  the  same  micropillar  array  as  (a).  Each  dark  fringe  corresponds  to  A/2 n  (A=405  nm  is  the 
wavelength  of  the  laser  and  n  =  1.3  is  the  refractive  index  of  water)  in  the  relative  thickness  of  the  liquid  film,  (c) 
Comparisons  of  the  simulated  and  measured  meniscus  along  the  red  line,  AB  with  the  ordinate  corresponding  to 
heights  relative  to  point  A. 
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where  AE  is  the  decrease  in  surface  energy  as  the  liquid  fills  one  unit  cell  and  A  V  is  the  corresponding 
volume  of  the  liquid  filling  one  unit  cell.  To  accurately  predict  the  change  in  surface  energy  and  volume, 
Surface  Evolver  (SE)[15]  was  used  to  simulate  the  shape  of  the  meniscus  with  minimum  surface  energy. 
A  unit  cell  consisting  of  four  quarter-pillars  filled  with  water  was  defined  as  shown  in  Fig.  3(a).  The 
effect  of  the  scalloped  features  were  included  in  the  simulations  as  an  effective  surface  energy  where  the 
scallops  were  approximated  as  a  series  of  semi-circles  with  a  roughness  factor,  rf=  7i/2[16].  The  intrinsic 
contact  angle  of  de-ionized  (DI)  water  on  the  silicon  surface  was  measured  by  a  goniometer  to  be  9  -  38°. 
Geometries  with  diameters,  d ,  and  periods,  /,  ranging  from  2  pm  to  8  pm  and  from  8  pm  to  15  pm, 
respectively,  were  considered. 

Fig.  3  compares  the  simulated  meniscus  shape  to  experimental  interference  measurements  on  a 
micropillar  array  with  diameter  d  -  6.4  pm  and  period  /  =  15  pm.  Fig.  3(a)  shows  the  three-dimensional 
meniscus  shape  obtained  by  SE  and  Fig.  3(b)  shows  the  top-down  view  of  the  experimentally  obtained 
interference  patterns  on  the  same  pillar  geometry.  Close  to  the  pillar  sidewalls  (<  1  pm),  accurate  data 
was  difficult  to  obtain  due  to  refraction.  Fig.  3(c)  compares  the  meniscus  shape  predicted  by  simulation 
and  obtained  experimentally  along  the  horizontal  (AB)  direction.  The  error  bars  in  the  data  were 
determined  based  on  the  width  of  the  dark  fringes.  The  experimental  results  were  well -predicted  by  the 
simulation.  Based  on  the  excellent  agreement  between  simulations  and  experiments,  SE  was  used  to 
obtain  general  correlations  for  the  volume  of  the  liquid,  AV,  and  the  area  of  the  menisci,  Am.  With  the 
prediction  of  meniscus  shape,  the  capillary  pressure  can  be  obtained  as 

_  AE  _  y-rfcos6-7rdh  +  ycos6(l2  -jrd2  / 4)- yAm 
cap  ~  AV  ~~  AV  ’ 

where  y  is  the  surface  tension  of  the  liquid,  8  is  the  intrinsic  contact  angle,  and  rf  is  the  roughness  factor 
on  the  sides  of  the  pillars.  The  first  and  second  term  represent  the  energy  change  as  the  liquid  wets  the 
sides  of  the  rough  pillar  and  the  smooth  bottom  surface  of  the  unit  cell,  respectively.  The  third  term 
corresponds  to  the  energy  change  due  to  the  formation  of  the  meniscus. 

However,  we  observed  in  our  experiments  that  two  distinct  time  scales  exist  in  the  propagation  process 
and  vary  as  functions  of  the  geometry.  The  time  required  for  the  liquid  front  to  sweep  from  one  pillar  to 
the  bottom  of  the  next  pillar  is  defined  as  rs  and  the  time  required  for  the  liquid  front  to  climb  and  wet  the 
side  of  a  pillar  is  defined  as  rp.  When  h>  /,  rs  and  rp  are  comparable,  such  that  the  averaged  approach 
described  above  to  determine  capillary  pressure  is  accurate.  On  the  other  hand,  when  h  <  /,  the  two  time 
scales  start  to  deviate  from  each  other  where  rs  becomes  dominant.  In  addition,  slight  variations  in  pillar 
diameter  and  local  wettability  become  important.  In  such  cases,  Eqn  (2)  is  no  longer  sufficient  to  describe 
the  driving  capillary  pressure.  We  plan  to  further  these  investigations  in  the  subsequent  year  to  obtain  a 
comprehensive  understanding  of  the  behavior  in  this  regime.  However,  these  sparse  pillar  geometries  are 
seldom  practical,  where  a  large  propagation  rate  is  often  desired. 

Subsequently,  we  determined  the  viscous  resistance  of  liquid  flow  in  micropillar  arrays  using 
Brinkmann’s  equation,  which  is  a  modified  form  of  the  Navier-Stokes  equation  with  a  permeability  term 
widely  used  in  porous  media  studies.  The  permeability  of  planar  pillar  arrays  was  obtained  using  the 
numerical  results  by  Sangani  and  Acrivos  for  dll  <  0.57  [17].  Conveniently,  pillar  arrays  with  dll  greater 
than  0.57  are  less  commonly  used  in  practice  due  to  the  large  viscous  losses  associated  with  such  densely 
packed  arrays.  Finite  element  simulations  were  carried  out  and  compared  to  our  model  and  show 
agreement  to  within  5%. 
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From  these  studies,  a  model  for  liquid  propagation  rates  was  developed,  and  verified  with  one¬ 
dimensional  propagation  experiments.  Comparisons  between  the  propagation  distance  as  a  function  of 
time  from  our  model  and  experiments  are  shown  in  Fig.  4(a)  for  d  =  2.9  pm,  /  =  5.5  pm,  and  h  -  8.3  pm, 
and  Fig.  4  (b)  for  d  =  2.82  pm,  /  =  8.1  pm,  and  h  =  8.3  pm.  Error  bars  were  obtained  by  repeating 
experiments  on  each  sample  at  least  three  times.  Experiments  were  also  performed  on  a  variety  of 
samples  with  d  =  1.8  -  5.3  pm,  l  =  4.5  -  15  pm,  and  h  =  8.3  -  25  pm.  Fig.  4  (c)  compares  the  propagation 
coefficients,  G,  from  experiments  and  the  model,  which  are  in  excellent  agreement  to  within  5%. 


d/I  Ratio 


(a)  (b) 

Fig.  5  Design  guidelines  for  the  optimization  of  pillar  array  geometries  for  h>l  and  dll  <  0.57.  (a)  Effective 
propagation  coefficients  predicted  by  the  model  for  varying  dll  and  hid  ratios,  (b)  Design  guidelines  for  maximum 
effective  propagation  coefficients  on  micropillar  arrays  with  various  aspect  ratios  and  the  corresponding  optimal  dll 
ratios. 


We  subsequently  developed  a  dimensionless  effective  propagation  rate  an  effective  propagation 
coefficient,  juG2  jyl,  was  obtained  that  is  functionally  dependent  on  only  the  dimensionless  geometrical 
parameters  hid  and  dll ,  where  /j  is  the  viscosity  of  the  liquid,  G  is  the  propagation  coefficient,  and  y  is  the 
surface  tension.  Fig.  5(a)  shows  the  relationship  for  water  on  silicon  pillars  with  various  dll  and  hid  ratios. 
The  relationship  is  valid  with  h  >  l  and  dll  <  0.57.  The  effective  propagation  coefficient  can  be  interpreted 
as  a  competition  between  the  viscous  resistance  and  the  driving  capillary  pressure.  As  the  ratio  of  hid 
increases,  the  effective  propagation  coefficient  increases  and  approaches  an  upper  theoretical  limit  as 
h — ►oo.  When  the  ratio  of  dll  increases,  the  viscous  resistance  and  the  capillary  pressure  both  increase  and 
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a  maximum  effective  propagation  coefficient  is  reached.  Beyond  this  limit,  the  viscous  resistance 
increases  faster  than  the  capillary  pressure,  leading  to  a  decrease  in  the  effective  propagation  coefficient. 
The  optimal  dll  ratio  and  the  corresponding  maximum  effective  propagation  coefficient  for  varying  hid 
ratios  are  shown  in  Fig.  5  (b).  The  outcomes  of  the  work  show  that  a  maximum  effective  propagation  rate 
can  be  achieved  by  optimizing  micropillar  geometries. 

2)  Evaporation-Induced  Non-Wetting  Droplets  on  Superhydrophilic  Surfaces 

We  discovered  during  our  experiments  that  liquid  droplets  could  exhibit  non-wetting  behavior  on 
superhydrophilic  microstructured  surfaces  at  superheated  conditions.  Such  a  non- wetting  state  is 
unfavorable  for  thin  film  evaporation  and  therefore,  we  investigated  the  detailed  mechanism  for  this 
behavior  and  outlined  the  conditions  to  avoid  non- wetting  of  liquids  in  thin  film  evaporation. 

We  studied  a  square  array  of  pillars  with  d  =  5.3  pm,  h  =  2l.l  pm,  /  =  27.5  pm,  roughness  ratio  rf=  1.48, 
and  solid  fraction  (j>s  -  0.03.  An  SEM  image  of  the  microstructured  surface  at  30°  inclination  is  shown  in 
Fig.  6a.  The  geometry  was  chosen  such  that  the  imbibition  condition  [18,  19]  was  met. 


(a) 

(b) 

10  nm 

5  mm 

Fig.  6  (a)  SEM  image  of  a  pillar  array  structure,  (b)  Image  of  front  and  back  side  of  a  typical  test  sample  used  in  the 
experiment. _ 


A  gold  resistive  heater  (Fig.  6b)  was  patterned  onto  the  backside  of  the  silicon  sample  using  ebeam 
evaporation  and  acetone  lift-off  processes.  A  DC  power  source  was  used  to  heat  the  sample  and  the 
temperature  was  measured  at  a  rate  of  10  Hz  using  a  T-type  thermocouple  (TC)  attached  to  the  backside 
of  the  test  sample  using  thermal  grease  (3G  Cool  Silver,  AI  Technology)  and  Kapton  tape.  Before  each 
experiment,  the  samples  were  cleaned  with  acetone,  methanol,  isopropanol  and  deionized  (DI)  water; 
dried  with  nitrogen  gas;  and  then  oxygen  plasma  treated  for  30  minutes  to  remove  organic  contaminants. 
After  the  plasma  cleaning  procedure,  the  surface  was  observed  to  be  superhydrophilic  and  droplets 
spontaneously  spread  upon  the  surface  with  near  zero  contact  angle.  Thoroughly  degassed  DI  water  was 
used  for  all  experiments  to  avoid  gas  bubble  nucleation.  A  pipette  was  used  to  gently  dispense  droplets 
with  volumes  between  15  and  25  pL  on  the  heated  substrate.  A  25  pm  diameter  gold  wire,  functionalized 
in  a  1  mM  ethanol  solution  of  1H,  1H,  2H,  2H-perfluorodecanethiol  (Sigma- Aldrich),  was  used  to  keep 
the  highly  mobile  droplets  in  position  for  visualization. 

The  superheat  required  to  deposit  a  non-wetting  droplet  on  top  of  the  microstructured  surface  was 
observed  to  be  higher  than  that  required  to  sustain  the  non-wetting  state  of  an  already  deposited  droplet. 
This  necessitated  a  strategy  that  would  allow  the  substrate  to  be  maintained  at  two  different  steady  state 
superheats:  a  higher  superheat  initially  and  a  relatively  lower  superheat  after  droplet  deposition.  This  was 
achieved  using  a  PID  controller  as  shown  in  the  schematic  circuit  setup  shown  in  Fig.  7  where  the 
substrate  was  initially  heated  directly  through  1 -2-3-4  to  establish  the  first  higher  steady  state  superheat 
and  through  1-5-6-7-3-4  to  establish  the  second  lower  steady  state  superheat.  The  second  steady  state 
superheat  is  achieved  by  toggling  a  3-way  switch  from  2  to  5  such  that  the  PID  controls  the  heating  power 
by  utilizing  a  closed  loop  temperature  feedback  system. 
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In  order  to  visualize  the  droplet  behavior  and  understand  the  mechanism  involved  during  wetting 
transition,  images  were  acquired  at  10  fps  which  are  later  synchronized  with  the  temperature  data  to 
extract  the  parameters  of  interest.  The  droplet  volume,  diameter  and  apparent  contact  angle  were 
calculated  by  processing  the  captured  images  through  a  custom  analysis  routine  written  in  MATLAB®.  To 
investigate  whether  or  not  a  water  droplet  wets  a  microstructured  surface  or  not  at  a  specified  superheat, 
ten  droplets  were  deposited  on  the  surface.  The  result  is  shown  on  Fig.  8  where  the  percentage  of  non¬ 
wetting  droplets  obtained  out  of  ten  experiments  is  shown  at  each  superheat.  All  droplets  deposited  below 
63  °C  superheat  wetted  the  surface  while  non- wetting  droplets  formed  at  all  times  when  the  superheat  was 
kept  higher  than  75  °C.  Between  63  °C  and  75  °C  superheat,  a  droplet  could  be  in  either  wetting  or  non¬ 
wetting  state  depending  on  how  gently  it  was  deposited  on  the  surface  [20].  However,  the  probability  of 
depositing  a  non- wetting  droplet  monotonically  increased  as  the  superheat  increased. 


A  typical  experiment  detailing  the  experimental  methodology  is  shown  in  Fig.  9.  A  steady  state  superheat 
of  94  °C  was  required  to  initially  deposit  a  non- wetting  droplet  while  only  26  °C  was  necessary  to  sustain 
the  non- wetting  state  of  an  already  deposited  droplet.  The  wetting  transition  occurred  at  26  °C,  35.1  s  after 
deposition. 
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(a) 


Fig.  9  (a)  Typical  time  lapse  images  of  an  evaporating  droplet  on  a  microstructured  surface  at  5  s  intervals,  (b) 
superheat  as  a  function  of  time,  (c)  corresponding  droplet  volume  as  a  function  of  time. 


To  explain  the  observations,  we  developed  a  simple  quasi-steady  1-D  lubrication  model  which  is  valid 
after  steady  state  superheat  was  achieved  (>18  s  on  Fig.  9b).  A  schematic  of  an  evaporating  droplet  on  top 
of  the  structured  surface  is  shown  on  Fig.  10.  The  temperature  at  the  droplet  base  (or  pillar  top)  is 
assumed  to  be  the  saturation  temperature  of  water  at  ambient  pressure.  Uniform  evaporation  throughout 
the  droplet  base  was  assumed.  The  apparent  contact  angle  of  the  droplet  was  estimated  using  the  Cassie- 
Baxter  equation  for  composite  surfaces  [21,  22],  a  porous  silicon  pillars  and  water  vapor  in  this  case, 


cos  6C  =  —1  +  0s(l  +  cos  Of), 


(3) 


where,  0f  and  6C  are  the  intrinsic  contact  angle  on  a  smooth  surface  and  the  apparent  contact  angle 
respectively 


We  further  assume  that  the  conducted  heat  is  completely  utilized  to  induce  uniform  evaporation  at  the 
liquid- vapor  interface  of  the  droplet  basal  area.  The  vapor  escapes  from  the  droplet  base  with  velocity  v  in 
the  z-direction  and  could  be  related  to  the  heat  conduction  through  the  composite  media  as 
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V  = 


(4) 


KffWy 

Pvap  hfg  h 


where  pvap  and  hfg  are  the  density  of  vapor  and  the  latent  heat  of  vaporization  respectively,  ATP  is  the 
temperature  difference  between  the  pillar  top  and  the  pillar  base,  and  kejf  is  the  effective  thermal 
conductivity  of  the  porous  media  composed  of  the  silicon  pillars  and  the  water  vapor  and  is  estimated  by 
Eq.  3  from  [23] 

^eff  —  (f  —  0s)^i?ap  "f  tps^sit  (5) 


where  kvap  and  ksi  are  the  thermal  conductivities  of  water  vapor  and  silicon  at  saturation  conditions  (ksi  at 
100  °C  ia  approximately  105  W/m.k)  [24]. 

Considering  that  h/R  «  1,  it  is  assumed  that  the  vapor  velocity  v  in  the  z-direction  is  small 
compared  to  the  radial  vapor  velocity  u.  The  continuity  equation  at  any  radial  position  under  the  droplet 
base  can  then  be  equated  with  Eq.  5  to  give 


keff^Tp  r 
^Pvaphfg  h 


(6) 


Further  assumptions  can  be  made  that  simplify  the  analysis  as  follows:  (1)  the  viscous  loss  is  dominated 
by  velocity  gradients  in  the  z-direction,  (2)  momentum  changes  in  the  axial  direction  can  be  neglected, 
and  (3)  the  pressure  depends  only  in  the  radial  direction.  The  above  assumptions  reduce  the  momentum 
equation  for  porous  media  to  Brinkman  equation  [25], 

d2u  s  dP  £U 

dz2  p  dr  +  K’  (7) 


where  s  is  the  porosity  and  K  is  the  permeability  of  the  porous  media  (unit  m2).  The  permeability  was 
calculated  using  an  asymptotic  expression  valid  for  flow  through  dilute  pillar  arrays,  (j)s  <  0.25  [26]. 
Equation  7  is  then  solved  by  applying  the  no-slip  boundary  conditions  at  the  pillar  base  and  pillar  top.  An 
average  radial  velocity,  which  depends  only  on  the  radial  direction,  can  then  be  found  by  integrating  the 
velocity  profile  from  the  pillar  base  to  the  tip.  This  average  radial  velocity  is  then  equated  with  the 
average  velocity  obtained  previously  from  Eq.  6  to  find  the  pressure  profile  along  the  radial  direction 
which  gives  rise  to  the  upward  anti- wetting  force  on  the  droplet  Fanti.wet  which  is  given  by 


HZPvap  ^eff  ^Tp 


r  anti-wet 


8 Pvaphfg 


2tanh 


h  17 

2^K 


\ 


hK 


© 


3/2 


/ 


R 4 
~h' 


(8) 


The  downward  wetting  force  on  the  droplet  Fwet  due  to  the  weight  of  the  droplet  and  the  surface  tension  is 
given  by 

F wet  PiiqVp  T  NndacosOf,  ^9) 

where  V  is  the  droplet  volume,  puq  is  the  density  of  liquid  water  at  saturation  temperature,  d  is  the  pillar 
diameter,  and  N  =  nR2ll 2  is  the  number  of  pillars  in  contact  with  the  droplet  base.  The  force  due  to  surface 
tension  is  the  dominant  wetting  force,  and  hence  the  scaling  Fwet  K R 2.  On  the  other  hand  the  anti-wetting 
force  scales  as  Fanti.wet  K  R4.  Hence,  due  to  evaporation,  the  anti-wetting  force  decreases  faster  than  the 
wetting  force.  The  superheat  required  to  suspend  a  droplet  on  top  of  the  pillar  array  can  be  obtained  by 
equating  the  two  competing  forces, 
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PuqVg  +  onNdcosOf 


A  Tp=- 


nzPvgp  feg// 

Bpvap  hfg 


2tanh 


Gr) 


3/2 


hK 

£ 


/ 


h 


(10) 


Equation  10  estimates  the  superheat  between  the  pillar  base  and  pillar  top  that  is  required  to  suspend  a 
droplet.  An  additional  temperature  drop  between  the  pillar  base  and  the  backside  of  the  substrate  is 
required  to  match  the  experimentally  measured  wall  temperature  with  the  model  prediction.  Incorporating 
this  additional  temperature  drop  across  the  sample  thickness  ( H  ~  620  pm),  the  overall  temperature  drop 
from  the  backside  of  the  sample  to  the  pillar  top  is 


AT  =  A  Tv 


Hkeff\ 
hksi  )' 


(11) 


Equation  11  provides  an  estimate  for  the  overall  superheat  AT  required  to  suspend  a  droplet  on  top  of  a 
microstructured  surface.  The  superheat  increases  as  the  droplet  decreases  in  size  ( ATaR '2),  consistent  with 
experimental  observations.  A  typical  experimental  data  and  the  model  prediction  are  shown  in  Fig.  11. 
The  simple  1-D  lubrication  type  model  qualitatively  captures  the  transition  mechanism  and  shows  a  good 
agreement  between  experimental  data  and  model  prediction. 


Fig.  1 1  Superheat  at  which  transition  happens  as  a  function  of  droplet  volume  predicted  by  our  model  and  validated 
by  experiments. _ 


In  summary,  we  demonstrated  that  a  superhydrophilic  surface  at  room  temperature  can  take  on  the 
characteristics  of  a  superhydrophobic  surface  if  the  substrate  is  superheated  beyond  a  critical  superheat. 
Experiments  performed  at  steady  state  conditions  on  superhydrophilic  surfaces  have  shown  that  a 
relatively  smaller  superheat  is  required  to  sustain  an  already  suspended  non-wetting  droplet  than  that 
required  when  initially  depositing  the  droplet.  We  have  also  demonstrated  experimentally  that  a  single 
surface  can  exhibits  both  extremes  of  wettability:  superhydrophilicity  and  superhydrophobicity.  The  1-D 
lubrication  type  model  developed  shows  good  agreement  with  experimental  data.  This  work  offers  new 
insights  into  the  design  of  structured  surfaces  for  phase-change  heat  transfer  applications  that  seek  high 
heat  flux  removal. 

3)  Pulsed  Evaporative  Transient  (PETT)  for  Temporally-resolved  Thermal  Property  Measurement 

Previous  studies  have  showed  that  the  heat  transfer  performance  of  thin  film  evaporation  is  highly 
dependent  on  the  microscale  liquid  film  morphology.  Numerical  simulations  have  compared  the  heat 
transfer  performance  of  different  microstructures  for  thin  film  evaporation  [27].  Meanwhile,  various 
experimental  approaches  have  been  used  to  study  the  relationship  between  meniscus  shape  and  heat 
transfer  behavior,  such  as  infrared  [28],  and  interference  imaging  [29].  However,  these  previous 
experimental  investigations  were  generally  steady  state  or  slowly  varying  measurements  where  a  constant 
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heat  flux  was  applied  and  a  steady-state  temperature  was  measured.  Accordingly,  the  average  heat 
transfer  coefficient  was  determined  over  a  long  period  of  time  (O(mins))  and  compared  with  imaging 
results.  While  these  approaches  are  of  practical  merit  where  the  actual  application  may  desire  steady-state 
performance,  to  fundamentally  understand  the  correlation  between  the  heat  transfer  coefficient  and  the 
liquid  film  morphology,  a  temporally-resolved  measurement  over  the  duration  of  the  evaporation  process 
is  of  interest. 

We  developed  a  new  metrology  technique,  Pulsed  Evaporative  Transient  Thermometry  (PETT),  to  obtain 
both  the  instantaneous  heat  capacity  and  heat  transfer  conductance,  i.e.,  the  product  of  the  heat  transfer 
coefficient  and  surface  area,  during  liquid  evaporation  from  micro  and  nanostructures.  This  approach 
experimentally  captures  the  transient  thermal  response  as  the  liquid  evaporates  from  micro  and 
nanostructures  by  a  series  of  applied  heat  pulses,  and  uses  a  lumped  system  model  to  extract  these  two 
quantities.  We  demonstrated  the  applicability  of  PETT  to  capture  the  evaporation  process  for  two  model 
systems,  copper  microwires  and  nanoporous  alumina  membranes  with  a  coated  metal  film.  In  the  case  of 
the  copper  microwires  (diameters~20  pm),  we  obtained  transient  heat  transfer  conductances  with  a 
temporal  resolution  of  1  second  and  captured  the  peak  heat  transfer  conductance  at  the  thin  film  stage  of 
1.3xl0"4  W/K,  i.e.,  -4*  as  high  as  the  dry  surface  (3.2xl0-5  W/K).  In  addition,  we  showed  that  the  heat 
transfer  behavior  were  distinct  between  the  oxidized  and  un-oxidized  surfaces,  which  correlated  well  with 
microscopic  visualizations  of  the  droplet  evaporation  process.  For  the  nanoporous  membranes,  a  higher 
temporal  resolution  of  0.2  seconds  was  achieved  owing  to  the  significantly  higher  heat  transfer 
conductance.  In  addition,  a  peak  in  the  heat  transfer  conductance  greater  than  1.2  W/K  during  the 
evaporation  process  was  recorded,  which  corresponds  to  the  thin  film  evaporation  state  with  the  smallest 
liquid  film  thickness.  This  PETT  technique  promises  to  be  a  valuable  quantitative  approach  to 
fundamentally  investigate  the  relationship  between  wetting  morphology  and  heat  transfer  behavior  for  a 
variety  of  micro  and  nanostructured  devices. 

The  experimental  setup  for  PETT  is  shown  in  Fig.  12a.  For  this  technique,  there  are  two  pre-requisites 
required  for  the  test  microdevice:  1)  The  Biot  number  of  the  device  ( Bi  =  hA/kL ,  where  hA  is  the  heat 
transfer  conductance,  k  is  the  thermal  conductivity,  and  L  is  the  characteristic  thickness  of  the  device) 
needs  to  be  low  (<0.1)  so  that  the  device  can  be  treated  as  a  lumped  system;  and  2)  The  device  needs  to 
have  a  finite  resistance  which  varies  with  temperature,  such  that  the  average  temperature  of  the  device  can 
be  measured  using  the  change  in  resistance.  This  test  device  is  connected  to  a  reference  resistor,  which  is 
stable  against  temperature  change,  and  a  function  generator,  which  provides  electrical  pulses.  After 
wetted  by  liquid,  the  microdevice  is  heated  by  applying  a  series  of  pulses  while  the  thermal  response  from 
the  resistance  change  is  captured  by  a  DAQ  system.  The  experimental  data  is  then  combined  with  the 
lumped  model  to  extract  the  heat  capacity  of  the  system  and  the  heat  transfer  conductance  described  as 
follows. 


When  the  device  is  heated  by  a  pulse,  the  temperature  response  is  determined  by 

^h  +  ^at  =  l  (i2) 

dt  c  c 

where  hA  is  the  heat  transfer  conductance  with  h  being  the  heat  transfer  coefficient  and  A  the  surface  area, 
C  is  the  heat  capacity  of  the  microdevice,  and  P  is  the  heating  power. 

By  solving  Eqn  (12),  the  temperature  difference,  AT,  as  a  function  of  time,  t ,  is  given  as 

AT  =  —  [1  -  e-^].  (13) 

hA  1  J 

In  the  case  of  microdevices  made  from  metal  or  coated  with  metal,  the  resistance  varies  linearly  with 
temperature.  Accordingly,  the  change  in  resistance  can  be  expressed  as 

aPRnr,  -—t-,  (14) 


AR  —  -——-2.  [1 

hA  L 


hA 

e^] 


where  a  is  the  temperature  coefficient  of  resistivity  and  R0  is  the  initial  resistance  of  the  device. 
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Therefore,  the  temporally-resolved  resistance  of  the  microdevice  is  determined  by  the  heat  transfer 
conductance,  M,  together  with  the  heat  capacity,  C.  Note  that  typically,  the  heat  transfer  conductance  of 
the  system  is  the  useful  quantity  to  determine  performance,  however  the  heat  transfer  coefficient  can  be 
determined  if  the  surface  area  of  the  microdevice  can  be  accurately  determined.  Fig.  12b  shows  two 
typical  responses  with  copper  micro  wires  (diameter  of  20  pm,  resistance  of  ~10  Q)  heated  by  pulses  with 
a  magnitude  of  600  mV  and  a  duration  of  0.2  seconds  with  water  droplets  evaporating  from  the  surface. 
As  the  droplets  evaporate,  the  heat  transfer  conductance  changes,  which  results  in  a  change  in  the 
resistance  response.  Meanwhile,  the  heat  capacity  of  the  system  can  be  obtained  using  the  response  at  the 

\lA 

early  stage  (/«—).  where  the  change  in  the  resistance  can  be  determined  as 

aPR0  hA  aPR0  (15) 

AR*-— - —  t  =  — - — t 

hA  C  C 

As  a  result,  the  change  in  the  resistance  is  approximately  linear  with  time  and  the  slope  indicates  the  heat 
capacity  as  shown  in  Fig.  12c. 


(a) 


Function 
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J  DAQcard  ]— 


Fig.  12  Pulsed  Evaporative  Transient  Thermometry  (PETT).  (a)  Schematic  of  the  experimental  setup.  A  function 
generator  provides  pulses  to  evaporate  the  liquid  on  the  surface  of  the  microdevice.  The  thermal  response  is 
recorded  as  a  function  of  time,  (b)  Example  changes  in  resistance  of  the  microdevice  with  time  during  a  heating 
pulse  of  0.2  seconds,  which  indicate  the  heat  transfer  conductance  (M).  (c)  Magnified  view  of  the  early  stage  of  the 
response  in  (b),  the  boxed  region,  where  the  heat  capacity  of  the  system,  C,  can  be  deduced  from  the  slope  of  the 
resistance  change  with  time.  These  data  were  obtained  from  the  microwire  (diameter  20  jam)  model  system  with 
water. 


Note  that  in  the  model,  the  heating  power  was  assumed  to  be  constant  to  achieve  the  ideal  exponential 
response.  In  actual  experiments,  a  constant  heating  power  is  difficult  to  implement.  Instead,  a  constant 
voltage  is  applied  to  the  reference  resistance  and  the  microdevice.  However,  by  matching  the  reference 
resistance  and  the  resistance  of  the  test  device,  the  variation  in  the  heating  power  can  be  minimized  as 
shown  below 


dP 

df 


d  (  u  y 

dT  \Rf  +  RJ  R 


dR  U2  2  R 

dT  (Rf  +  R)2  ^  Rf  +  ^ 


(16) 


where  Rf  is  the  reference  resistance  and  U  is  the  total  voltage.  Therefore,  by  choosing  Rf  ~  R ,  the  variation 
in  the  heating  power  can  be  neglected  and  the  proposed  model  will  accurately  capture  the  experiments. 

We  first  used  copper  micro  wires  as  a  model  system  to  demonstrate  the  capability  of  the  PETT  technique, 
where  we  examined  the  difference  in  evaporation  behavior  and  corresponding  heat  transfer  conductances 
on  oxidized  and  un-oxidized  microwires.  Copper  wires  were  procured  (McMaster  Carr)  with  diameters 
of  76  pm,  and  were  subsequently  etched  by  iron(III)  chloride  solution,  decreasing  the  diameter  to  ~20 
pm,  and  increasing  the  electrical  resistance  of  the  sample  to  the  range  of  10-30  Q.  With  these  resistance 
values,  a  simple  2-wire  configuration  could  be  used  to  achieve  accurate  resistance  measurements.  After 
etching,  we  oxidized  some  of  the  copper  microwires  to  effectively  increase  the  hydrophilicity  of  the 
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copper  surface  [30].  The  scanning  electron  microscope  (SEM)  images  of  the  oxidized  surfaces  are  shown 
in  Fig.  13(a).  The  oxidation  process  significantly  enhanced  both  the  micro  and  nanoscale  roughness  of  the 
micro  wire  in  comparison  with  the  un-oxidized  micro  wires  in  Fig.  13(b). 


The  microwires  were  placed  into  the  experimental  setup  in  Fig.  12.  PETT  measurements  were  first 
obtained  on  dry  micro  wires  to  validate  the  method  and  determine  the  uncertainty  of  the  results.  The 
temporally  resolved  heat  capacity  and  heat  transfer  conductance  on  dry  oxidized  and  un-oxidized  copper 
wires  are  shown  in  Fig.  13(c)  &  (d),  respectively.  On  oxidized  micro  wires,  the  average  heat  capacity  and 
heat  transfer  conductance  over  50  seconds  were  approximately  3.7xl0"5  J/K  and  2.1xl0"4  W/K,  with 
uncertainties  of  ±1.7xl0"6  J/K  and  ±6.3xl0"6  W/K,  respectively.  On  un-oxidized  microwires,  the  average 
heat  capacity  and  heat  transfer  conductance  over  50  seconds  were  approximately  7.5xl0"5  J/K  and 
4.3xl0"4  W/K,  with  uncertainties  of  ±1.6xl0"5  J/K  and  ±1.3xl0"5  W/K,  respectively.  The  uncertainties 
were  one  order  of  magnitude  smaller  compared  to  the  average  values,  indicating  good  consistency  and 
reliability.  Note  that  the  measured  heat  capacity  and  heat  transfer  conductance  on  the  un-oxidized  wires 
were  ~2*  of  that  on  the  oxidized  wire.  This  was  due  to  the  fact  that  the  oxidized  wires  used  in  the  tests 
were  approximately  half  the  length  of  the  un-oxidized  ones.  We  chose  the  different  lengths  in  experiments 
for  practical  reasons.  Ideally,  the  wires  should  be  as  long  as  possible  so  that  the  electrical  resistance  can 
be  large  for  easier  measurements.  However,  during  the  oxidation  process,  the  samples  were  placed  in  the 
oxidation  solution  with  agitation.  Long  wires  with  such  small  diameters  (~20  pm)  tended  to  break  more 
easily  due  to  agitation.  Therefore,  we  reduced  the  lengths  of  the  wires  for  oxidation. 

Afterwards,  we  wetted  the  surface  with  water  and  captured  the  evaporation  process  without  heating  using 
a  CMOS  camera  (Phantom  V7.1,  Vision  Research  Inc.)  at  a  frame  rate  of  10  frames  per  second  through  a 
white  light  microscope  (Eclipse  LV-100,  Nikon)  with  a  20x  objective  (NA  =  0.60). 
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Fig.  14  PETT  measurements  during  evaporation  on  an  oxidized  copper  microwire,  (a)  White-light  microscope 
images  of  a  typical  evaporation  process  of  a  liquid  water  droplet  on  an  oxidized  wire.  Strong  pinning  of  the  droplet 
resulted  in  a  thin  liquid  film  covering  the  surface  at  the  late  stage  of  evaporation  (t=  12.5  s,  bright  reflection 
highlighted  by  the  circle),  (b)  Temporally-resolved  heat  capacity  (C)  and  heat  transfer  conductance  (M)  on  the 
wetted  wire  with  a  resolution  of  1  second.  A  peak  in  the  heat  transfer  conductance  corresponds  to  the  high  heat 
transfer  rate  associated  with  the  thin  film  evaporation  stage. 


The  time-lapse  images  in  Fig.  14(a)  show  that  the  liquid  droplet  was  pinned  due  to  the  enhanced 
roughness  and  hydrophilicity  and  the  thickness  of  the  droplet  decreased  during  the  evaporation  process. 
Therefore,  at  the  late  stage  of  the  evaporation  (t=  12.5  s  in  this  specific  experiment),  thin  liquid  films  were 
formed,  as  highlighted  by  the  circle  in4.  18(a).  Based  on  these  visualizations,  a  peak  in  heat  transfer 
conductance  was  expected  corresponding  to  this  thin  film  stage.  The  PETT  measurement  during  an 
example  evaporation  experiment  is  shown  in  Fig.  14(b).  Each  heating  pulse  lasted  for  0.2  seconds  to 
achieve  a  significant  increase  in  resistance.  A  subsequent  0.8  seconds  waiting  period  between  pulses  was 
used  in  order  for  the  sample  to  cool  to  the  ambient  temperature  prior  to  the  next  pulse.  As  a  result,  the 
overall  temporal  resolution  was  1  second.  The  heat  capacity,  C,  of  the  device  decreased  from  the  initial 
value  of  -l.OxlO"4  J/K  as  liquid  evaporated  to  a  steady  dry  value  of  -5x1  O'5  J/K  at  the  end  of  the 
experiment.  The  heat  transfer  conductance,  M,  increased  to  reach  a  peak  value  of  -1.3  xlO"3  W/K  and 
then  decreased  to  the  dry  value  of  -3.2xl0"4  W/K  towards  the  end  of  the  experiment,  as  expected. 
Accordingly,  the  peak  heat  transfer  conductance  was  4x  as  large  as  that  when  the  wire  was  dry. 

For  comparison,  experiments  were  also  carried  out  on  un-oxidized  copper  microwires  with  smooth 
surfaces.  Time-lapse  images  of  a  representative  evaporation  process  on  the  un-oxidized  wires  are  shown 
in  Fig.  15(a).  Due  to  the  reduced  roughness  and  consequently  less  pinning  of  the  droplets,  the  base  area  of 
droplets  decreased  continuously  during  the  evaporation  process  and  the  thin  film  stage  was  not  observed. 
Therefore,  no  peak  in  heat  transfer  conductance  was  expected,  as  confirmed  by  the  temporally -resolved 
measurement  of  heat  capacity  and  heat  transfer  conductance  in  Fig.  15(b).  Similar  with  the  oxidized 
wires,  the  heat  capacity  gradually  decreased  from  -1.5  xlO"4  J/K  to  -7xl0"5  J/K  as  water  evaporated. 
However,  the  heat  transfer  conductance  also  decreased  gradually  from  -1.5xl0"3  W/K  to  -0.5xl0"3  W/K 
with  no  peak.  The  heat  transfer  behavior  was  distinctly  different  between  the  evaporation  processes  on 
oxidized  and  un-oxidized  copper  wires  due  to  the  difference  in  surface  property  and  liquid  morphology. 
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Fig.  15  PETT  measurements  during  evaporation  on  an  un-oxidized  copper  microwire,  (a)  White-light  microscope 
images  of  a  typical  evaporation  process  of  a  liquid  water  droplet  on  an  un-oxidized  microwire.  The  wetted  area 
shrunk  gradually  during  the  evaporation  process,  (b)  Temporally-resolved  heat  capacity  (C)  and  heat  transfer 
conductance  (M)  on  the  wetted  wire  with  a  resolution  of  1  second.  There  is  no  peak  in  the  heat  transfer  conductance 
because  the  high  heat  transfer  rate  with  the  thin  film  evaporation  stage  was  not  observed. 


Using  microwires,  we  showed  that  the  PETT  technique  is  capable  of  measuring  transient  thermal 
properties.  While  the  microwires  offer  a  convenient  platform  for  fundamental  investigations,  they  are  not 
common  configurations  for  practical  applications.  Moreover,  the  low  heat  transfer  conductance  requires 
long  cooling  period  between  pulses,  which  limits  the  temporal  resolution.  We  further  expanded  the  PETT 
technique  to  a  more  practical  configuration  for  thin  film  evaporation  using  nanoporous  membranes  wetted 
by  isopropyl  alcohol.  Isopropyl  alcohol  was  used  in  this  case  to  allow  easy  wetting  of  the  membrane  due 
to  the  low  surface  tension,  and  to  delay  electrolysis  when  in  contact  with  the  metal  electrode.  The  large 
surface  area  and  complete  wetting  allowed  significantly  higher  heat  transfer  conductances,  which  greatly 
reduce  the  required  cooling  period.  Therefore,  we  targeted  0.2  seconds  of  temporal  resolution  with  a  pulse 
duration  of  0.1  seconds  and  cooling  period  of  0.1  seconds.  As  will  be  discussed  in  more  detail  later, 
however,  there  is  a  trade-off  between  obtaining  high  temporal  resolution  and  high  signal-to-noise  ratio. 

The  average  pore  diameter  of  the  procured  alumina  membranes  (Synkera  Inc.)  was  150  nm  with  a 
thickness  of  100  pm.  A  30  nm  thick  layer  of  platinum  was  sputtered  onto  the  membrane  surface,  which 
served  as  both  the  heater  and  an  area-averaged  temperature  sensor.  Images  of  the  coated  membrane  are 
shown  in  Fig.  16(a).  Typical  changes  in  resistance  with  time  responses  are  shown  in  Fig.  16  (b)  for 
different  heat  transfer  conductances.  Due  to  the  large  heat  transfer  conductance  and  short  pulse 
introduced,  however,  the  increase  in  resistance  during  the  pulse  (<0.04  Q)  became  much  smaller 
compared  to  that  in  the  microwire  configuration  (>0.2  Q).  As  a  result,  noise  from  data  acquisition  system 
became  much  more  significant  compared  to  the  desired  signal. 

The  temporally-resolved  heat  transfer  conductance  measured  by  PETT  is  shown  in  Fig.  16(c).  Initially 
when  the  membrane  was  completely  wetted  and  flooded  with  IPA,  the  heat  transfer  conductance  was 
steady  around  0.38  W/K.  As  evaporation  occurred,  a  clear  peak  in  heat  transfer  conductance  was  observed 
prior  to  dry-out,  which  corresponds  to  the  stage  where  the  membrane  was  wetted  by  a  thin  liquid  film. 
Later,  as  the  thin  liquid  film  evaporated,  the  heat  transfer  coefficient  decreased  to  0.14  W/K,  which  is 
associated  with  natural  convection  and  conduction  loss  through  the  electrical  contacts. 

The  changes  in  resistance  during  the  peak  in  heat  transfer  conductance  are  shown  as  the  red  data  points  in 
Fig.  16(b).  A  curve  fit  to  the  data  yielded  a  heat  transfer  conductance  of  3.6  W/K  with  an  increase  of 
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0.006  Cl  in  the  resistance,  which  is  smaller  than  the  standard  deviation  of  0.01  Cl.  This  fitted  result 
corresponds  to  an  upper  bound  of  heat  transfer  conductance,  plotted  in  Fig.  16c  as  a  red  square.  It  is 
possible  that  the  actual  increase  in  the  resistance  was  higher  but  was  masked  by  the  noise.  However,  the 
upper  bound  for  the  actual  increase  should  be  twice  that  of  the  standard  deviation  (<0.02  Q).  Based  on 
Eqn  (3),  the  upper  bound  of  the  actual  resistance  increase  corresponds  to  a  lower  bound  in  heat  transfer 


conductance,  which  is  ( [hA)min  =  rj*P.R°  =  1.2  W/K.  This  value  is  also  plotted  in  Fig.  16(c)  as  a  red 

(A  R)max 

dot.  The  lower  bound  in  heat  transfer  conductance  at  the  peak  remains  one  order  of  magnitude  larger  than 
that  at  the  dry  stage.  Despite  the  uncertainty  in  the  actual  peak  heat  transfer  conductance,  the  results 
indicate  that  the  peak  in  the  thin  film  evaporation  stage  exists.  Due  to  the  significant  noise  and  small 
resistance  change,  however,  we  were  not  able  to  determine  the  heat  capacity. 


(b) 


Fig.  16  PETT  measurements  on  a  nanoporous  membrane,  (a)  Images  of  the  nanoporous  membrane  (left)  with  a  30 
nm  thick  platinum  layer  deposited  on  the  membrane.  SEM  image  (right)  shows  the  nanoscale  pore  geometries  after 
platinum  deposition,  (b)  Example  responses  of  the  change  in  resistance  with  time  during  evaporation  with  low  and 
high  heat  transfer  conductance.  The  change  in  resistance  was  much  smaller  than  that  in  the  microwire  system  due  to 
the  shorter  (0.1  seconds)  pulses  and  larger  heat  transfer  conductance  (>0.1  W/K).  Consequently,  the  lower  signal- 
to-noise  ratio  led  to  a  much  higher  measurement  uncertainty.  A  clear  distinction  in  the  responses,  however,  was  still 
observed  between  the  two  cases  with  different  heat  transfer  conductances,  (c)  Temporally-resolved  heat  transfer 
conductance.  The  peak  corresponds  to  the  thin  film  evaporation  stage.  An  upper  bound  of  3.6  W/K  (red  square)  at 
the  peak  was  determined  by  curve  fitting  and  a  lower  bound  of  1.2  W/K  (red  dot)  was  estimated  using  the  maximum 
temperature  rise  during  the  pulse. 


Based  on  the  results  from  the  two  types  of  model  systems,  we  show  that  there  are  trade-offs  in 
implementing  this  technique  associated  with  temporal  resolution  and  accuracy.  We  showed  that  for 
systems  with  high  heat  transfer  conductance  (membrane),  the  cooling  period  between  the  applied  pulses 
can  be  much  shorter  to  achieve  improved  temporal  resolution,  but  the  uncertainty  of  the  results  can 
become  significant.  Meanwhile,  for  systems  with  low  heat  transfer  conductance  (microwires),  the  large 
increase  in  resistance  allows  for  higher  signal-to-noise  ratios  and  more  accurate  measurements  of 
transient  heat  capacity  and  heat  transfer  conductance;  however,  the  temporal  resolution  is  limited. 
Accordingly,  the  considerations  of  accuracy  and  temporal  resolution  need  to  be  thoroughly  considered 
before  determining  the  proper  parameters  to  implement  the  PETT  technique.  In  addition,  improving  the 
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signal  to  noise  ratio  in  the  DAQ  system  would  also  be  helpful  to  increase  the  accuracy  of  the 
measurements. 

In  summary,  a  new  metrology  technique,  Pulsed  Evaporative  Transient  Thermometry  (PETT),  was 
presented  to  determine  instantaneous  heat  capacity  and  heat  transfer  conductance  of  the  system.  This 
technique  utilizes  short  pulses  to  heat  the  sample  and  captures  the  transient  response  to  extract  these 
quantities.  We  demonstrated  this  approach  to  two  model  systems  of  copper  micro  wires  with  different 
surface  treatments  and  nanoporous  membranes.  With  different  system  properties,  optimized  accuracy  and 
temporal  resolution  can  be  achieved  by  tuning  the  duration  of  heating  pulses  and  the  cooling  period.  This 
approach  promises  wide  application  to  a  variety  of  micro  and  nanostructures  to  improve  understanding  of 
phase  change  heat  transfer  processes. 

4)  Thin  Film  Evaporation  from  Nanoporous  Membranes 

In  addition  to  using  microstructured  surfaces,  we  realized  upon  the  course  of  this  study  that  a  nanoporous 
membrane  configuration  promises  the  highest  performance.  Accrdoingly,  we  performed  thin  film 
evaporation  experiment  from  nanoporous  membranes  using  isopropyl  alcohol  (IP A)  as  the  coolant.  With 
average  pore  diameters  of  150  nm,  absolute  negative  pressures  as  low  as  ~-300  kPa  were  demonstrated. 
These  pressures  enabled  interfacial  heat  fluxes  of  ~96  W/cm2  using  IPA  as  the  working  fluid. 
Furthermore,  we  provide  design  guidelines  to  achieve  heat  fluxes  over  1000  W/cm2  by  creating 
nanoporous  membranes  with  thicknesses  on  the  order  of  1  pm.  This  study  is  a  first  step  towards 
demonstrating  the  feasibility  of  a  high  performance  heat  flux  device  using  thin  film  evaporation. 

The  nanoporous  membranes  used  in  our  experiment  was  anodized  aluminum  oxide  membranes  (Unikera, 
Synkera  Inc.)  with  diameter  of  13  mm.  The  average  pore  diameter  was  150  nm.  A  35  nm  thick  platinum 
layer  was  deposited  on  top  of  the  membrane  by  sputtering,  which  was  used  as  resistance  heater  as  well  as 
temperature  sensor.  The  liquid  was  driven  by  capillary  pressure  as  well  as  disjoining  pressure  through  the 
pores  to  form  the  menisci.  An  example  macroscale  image  and  SEM  image  of  the  membrane  after  metal 
deposition  are  shown  in  Fig.  17(a)  and  (b),  respectively.  Note  that  the  average  pore  diameter  of  150  nm 
and  porosity  of  32%  is  based  on  manufacturer’s  specifications.  Based  on  image  analysis  of  the  SEM 
image  after  metal  deposition,  the  porosity  of  this  membrane  is  reduced  to  14%  by  the  metal  layer.  Later, 
we  show  that  the  porosity  after  metal  deposition  ranges  and  needs  to  be  quantified  for  every  batch  of 
fabricated  membranes. 


(a)  (b)  (c) 


Fig.  17  (a)  Macroscale  image  of  the  membrane  after  metal  deposition,  (b)  and  (c)  SEM  image  of  the  membrane  after 
metal  deposition.  The  porosity  for  this  specific  membrane  was  reduced  to  14%. 

During  the  experiment,  the  temperature  was  measured  by  the  change  in  the  resistance  of  the  metal  layer. 
However,  due  to  the  nanoscale  geometry  of  the  platinum  layer  and  the  sputtering  fabrication  process,  the 
electrical  property  of  the  metal  layer  is  different  from  that  of  the  bulk  material.  Therefore,  the  temperature 
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coefficient  of  resistivity  was  calibrated  by  an  oven  and  a  RTD  thermal  probe.  The  result  of  the  calibration 
is  shown  in  Fig.  18.  The  temperature  coefficient  of  resistance  was  found  to  be  0.00143  K"1. 
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Fig.  18  Calibration  of  the  temperature  coefficient  of  resistivity  for  the  sputtered  metal  layer.  The  temperature 
coefficient  was  found  to  be  0.00143  K"1. 


In  the  experiment,  the  membrane  was  attached  to  a  liquid  feed  column  with  the  edge  of  the  membrane 
sealed  by  vacuum  epoxy  (Torr  Seal,  Varian  Inc.)  to  ensure  that  liquid  can  only  transport  through  the 
nanopores  in  the  membrane,  as  shown  in  Fig.  19(a).  The  area  that  can  be  wetted  by  the  fluid  in  the 
membrane  is  0.5  cm2.  Considering  the  porosity  of  14%,  the  actual  interfacial  area  available  for 
evaporation  is  0.07  cm2.  The  liquid  feed  column  is  placed  in  a  vacuum  chamber  (Kurt  J.  Lesker).  A 
vacuum  pump  removes  the  vapor  generated  during  evaporation  to  eliminate  the  thermal  resistance 
associated  with  vapor  transport.  The  pressure  of  the  chamber  is  monitored  by  a  pressure  transducer.  The 
bottom  of  the  liquid  feed  column  is  connected  to  a  liquid  reservoir  outside  the  vacuum  chamber  at 
atmosphere.  A  voltage  was  supplied  to  the  metal  layer  to  provide  the  heat  fluxes.  The  voltage  and  current 
was  measured  by  multimeters  (Multimeter  2001,  Keithley)  and  recorded  by  a  computer  using  Lab  View 
(National  Instruments)  program  such  that  the  total  heating  power  and  temperature  can  be  recorded.  The 
data  was  acquired  at  a  frequency  of  1  Hz.  The  schematic  and  actual  image  of  the  setup  is  shown  in  Fig. 
19(b)  and  (c),  respectively.  At  each  heating  power,  the  membrane  is  allowed  to  reach  steady  state  where 
the  temperature  variation  is  typically  lower  than  2  degrees  over  100  seconds.  After  that,  the  heating  power 
is  varied  to  gather  more  data  points. 


(a)  (b)  (c) 

Fig.  19  Experimental  setup  (a)  Image  of  the  membrane  attached  to  the  liquid  feed  column,  (b)  Schematic  showing 
the  experimental  setup,  (c)  Image  of  the  experimental  setup  placed  in  the  vacuum  chamber. 
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In  the  experiments,  the  heat  was  dissipated  by  evaporation,  conduction  through  the  fixture,  and 
convection  through  the  liquid.  The  conduction  heat  loss  through  the  fixture  can  be  determined  by 
measuring  the  temperature  change  as  a  function  of  heating  power  without  liquid  supply,  as  shown  in  Fig. 
20.  In  the  test,  the  chamber  was  pumped  down  the  actual  experiment  condition  (~60  Pa).  The  heat  loss 
through  the  fixture  was  found  to  be  around  0.01  W/K. 


Fig.  20  Temperature  of  the  membrane  as  a  function  of  heating  power  without  liquid  supply.  The  conduction  loss 
through  the  fixture  was  found  to  be  around  0.01  W/K. _ 


The  heat  loss  through  convection  through  the  liquid  is  estimated  using  the  correlation  given  by  Kadambi 
and  Drake  for  external  natural  convection  of  a  heated  horizontal  plate  facing  down  [31] 


Nut  =  0.82  Ra 


1/5 


(17) 


lL  ~ 

where  RaL  is  the  Rayleigh  number  determined  by 

RaL=—  ATL3  (18) 

va 

Given  the  fluid  properties  of  isopropyl  alcohol,  the  heat  loss  due  to  convection  can  be  estimated  to  be 

QtA)  =  0.006  W/K  (19) 


Note  that  the  correlation  used  here  is  for  external  natural  convection  where  the  fluid  is  not  confined 
compared  to  the  practical  experiment  fixture  where  the  fluid  is  confined  in  the  feed  column.  Therefore  the 
calculated  heat  loss  should  be  an  overestimation  compared  to  the  actual  case. 

The  typical  temperature  rise  during  the  experiment  is  around  30-50  K.  Therefore,  the  maximum  total  heat 
loss  will  be  <0.8  W,  which  is  very  small  compared  to  the  applied  heating  power. 


The  variation  of  the  steady  state  temperature  of  the  membrane  with  the  increase  in  heat  flux  is  shown  in 
Fig.  21(a).  Different  symbols  indicate  different  tests,  where  good  consistency  was  observed.  Data  was 
acquired  using  isopropyl  alcohol  (IP A)  as  the  coolant.  We  used  IPA  instead  of  water  due  to  the  low 
surface  tension  of  IPA,  which  allows  wetting  of  the  membrane  even  when  the  membrane  gets 
contaminated  and  become  less  hydrophilic.  Future  effort  will  be  focused  on  fabricating  surfaces  with 
robust  hydrophilicity.  As  shown  in  Fig.  21(a),  the  initial  heat  transfer  coefficient  was  relatively  low  with 
low  heat  flux.  This  result  is  because  at  the  early  stage,  the  liquid  flow  driven  by  the  pressure  difference 
applied  by  the  vacuum  pumping  is  higher  than  the  evaporation  rate.  The  top  surface  of  the  membrane  was 
flooded  by  the  excessive  amount  of  liquid,  which  introduced  a  large  thermal  resistance.  As  the  heat  flux 
was  increased,  the  evaporation  rate  became  higher  than  the  flow  rate  driven  by  pressure  difference  and  the 
excessive  amount  of  liquid  on  the  surface  on  the  membrane  evaporated.  The  heat  transfer  coefficient 
increased  significantly  as  the  liquid  menisci  receded  back  into  the  pores  and  formed  thin  film  regions. 
With  further  increase  in  the  heat  flux,  the  menisci  shape  yielded  larger  thin  film  regions  resulting  in  the 
continuous  increase  of  the  heat  transfer  coefficient.  At  this  stage  with  higher  heat  transfer  coefficient,  the 
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liquid  was  driven  through  the  pores  by  a  combination  of  the  vacuum  pumping  pressure,  the  capillary 
pressure  and  the  disjoining  pressure.  However,  with  further  increase  in  the  heat  flux,  the  total  driving 
pressure  was  no  longer  able  to  supply  enough  liquid  for  evaporation,  and  the  membrane  started  to  dry  out 
with  significant  rise  of  the  membrane  temperature.  The  images  showing  the  flooding  and  the  thin  film 
evaporation  states  of  the  membrane  are  shown  in  Fig.  21(b)  and  (c),  respectively.  The  heat  transfer 
coefficient  of  thin  film  evaporation  can  be  determined  as  1.3  W/cm2K  using  the  slope  of  the  data  points  in 
the  thin  film  operating  state  (dashed  line  in  Fig.  21(a)).  Before  the  drying  out  occurs,  the  maximum 
interfacial  heat  flux  was  -50  W/cm2  based  on  the  actual  interfacial  area  of  0.07  cm2  at  the  working 
temperature  of  345  K. 
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Fig.  21  Experimental  results  on  100  pm  membranes,  (a)  Relation  between  total  dissipated  power  and  temperature.  A 
drop  in  the  temperature  was  observed  at  the  transition  between  “flooding”  and  thin  film  regime.  The  heat  transfer 
coefficient  was  determined  to  be  1.3  W/cm2K  and  the  maximum  interfacial  heat  flux  was  50  W/cm2  based  on  the 
actual  pore  area  of  0.07  cm2,  (b)  Image  of  “flooding”  of  the  membrane.  Large  droplets  were  observed  on  top  of  the 
membrane  which  limits  the  heat  transfer  coefficient,  (c)  Image  of  the  thin  film  regime  of  the  membrane.  The 
droplets  disappeared. 


The  transition  between  flooding  and  thin  film  evaporation  occurred  at  a  heating  power  of  approximately 
1.2-1. 3  W,  where  the  fluid  transport  through  the  nanopores  is  purely  driven  by  the  pressure  difference 
(atmospheric  pressure,  -lOOkPa)  across  the  membrane.  While  at  the  dry-out  point  (-3.9  W),  the  fluid  is 
driven  by  the  combination  of  atmospheric  pressure  and  capillary  pressure.  Based  on  energy  conservation, 
the  flow  rate  before  the  dry-out  should  be  3. 9/1. 3=3  times  as  large  compared  to  the  flow  rate  at  the 
transition  point  between  flooding  and  thin  film  evaporation.  Therefore,  the  total  driving  pressure 
including  atmospheric  pressure  and  capillary  pressure  should  be  -300  kPa.  As  a  result,  the  fluid  near  the 
meniscus  is  in  negative  pressure  (—200  kPa).  This  negative  pressure  is  allowed  by  nanoscale  confinement 
of  the  pores  which  prevents  the  cavitation  of  fluid  [32].  Using  Young-Laplace  equation,  the  maximum 
capillary  pressure  can  be  estimated  as 
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(20) 


2  a 

PCap  =  —  =  746  kPa 

where  o  is  the  surface  tension  of  the  liquid  (-0.019  at  345K)  and  r  is  the  average  radius  of  the  pores  (-50 
nm  after  metal  deposition). 

The  actual  demonstrated  capillary  pressure  is  significantly  lower  than  the  theoretical  limit  due  to  a  few 
reasons.  First  of  all,  Young-Laplace  equation  assumes  a  circular  pore  while  the  actual  pores  in  the 
membrane  are  not  circular  and  have  a  distribution  in  the  pore  size,  as  can  be  seen  in  the  SEM.  Secondly, 
the  fluid  used  in  the  experiment  is  only  99%  pure.  Significant  contamination  was  observed  after  the 
experiment,  as  shown  in  Fig.  22(a)  and  (b).  Such  contamination  reduces  the  capillary  pressure  available  to 
drive  the  liquid  transport.  _ 


(a)  (b) 

Fig.  22  Contamination  of  membranes  during  evaporation,  (a)  Image  of  the  membrane  before  experiments,  (b)  Image 
of  membranes  after  the  experiments.  The  arrows  highlight  the  contaminated  marks. 


In  order  to  evaluate  the  demonstrated  interfacial  heat  flux  of  50  W/cm2,  the  critical  heat  flux  (CHF)  of 
IPA  in  pool  boiling  was  estimated  using  Kandlikar’s  model  for  a  benchmark  comparison  [33].  The  CHF 
is  dependent  on  the  orientation  of  the  heated  surface  and  ranges  from  39  W/cm2  (vertical  surface)  to 
73  W/cm2  (horizontal  surface).  The  demonstrated  heat  flux  is  comparable  with  CHF  of  pool  boiling.  The 
major  limitation  in  the  heat  flux  is  the  supply  of  fluid.  Due  to  the  relatively  large  thickness  (100  pm)  of 
the  membrane,  the  viscous  resistance  of  liquid  flow  through  the  nanopores  is  significant.  However, 
assuming  a  constant  driving  capillary  pressure  of  200  kPa  and  same  pore  geometries,  the  maximum  heat 
flux  would  be  inversely  proportional  to  the  thickness  of  the  membrane.  Therefore,  the  decoupling 
between  capillarity  and  viscous  resistance  and  significantly  higher  heat  fluxes  can  be  achieved  by 
reducing  the  thickness  of  the  membrane. 

In  order  to  validate  the  scaling,  experiments  using  thinner  membranes  were  performed.  The  thickness  of 
the  new  membrane  is  reduced  to  50  pm.  SEM  image  of  the  pore  geometries  is  shown  in  Fig.  23  (a),  by 
which  the  porosity  after  the  deposition  was  determined  as  24%  and  the  total  pore  area  is  0.12  cm2.  The 
results  are  summarized  in  Fig.  23(b).  The  pressure  in  the  chamber  during  the  experiment  is  around  200- 
300  Pa.  Due  to  the  reduced  thickness  and  higher  porosity,  a  much  larger  amount  of  fluid  can  be  pushed 
through  the  membrane  at  low  heating  power  to  form  a  liquid  pool  on  top  of  the  membrane  where 
nucleation  was  observed  in  contrast  to  the  liquid  droplets  observed  on  the  100  pm  thick  membranes.  As  a 
result,  significant  temperature  fluctuation  was  seen  at  low  heating  power  which  prohibits  recording 
accurate  data.  Meanwhile,  with  increased  heating  power,  the  flooding  can  be  removed  and  the 
temperature  becomes  stabilized.  The  comparison  between  low  heating  power  and  high  heating  power  is 
shown  in  Fig.  23(b).  Therefore,  in  these  experiments,  the  data  points  at  low  heating  power  are  not  shown. 
The  transition  from  flooding  and  thin  film  evaporation  is  determined  as  the  point  where  the  fluctuation  in 
temperature  starts  to  diminish,  which  is  around  2.9  W. 
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Fig.  23  Experimental  results  on  50  pm  thick  membranes,  (a)  SEM  images  of  platinum  coated  nanoprous  membranes 
with  porosity  of  24%  and  thickness  of  50  pm.  (b)  Relation  between  dissipated  power  and  temperature.  The  heat 
transfer  coefficient  was  determined  to  be  1.8  W/cm2K  and  the  maximum  interfacial  heat  flux  was  96  W/cm2.  (c) 
Typical  area-averaged  temperature  responses  of  flooding  and  thin  film  evaporation  regimes  over  100  seconds.  The 
temperature  fluctuations  were  large  during  flooding  (-15  K)  with  2.3  W  heating  (bottom)  and  small  (-1  K)  during 
thin  film  evaporation  with  -3.6  W  heating  (top) 


Using  the  slope  of  the  data  points  in  the  thin  film  region,  the  heat  transfer  coefficient  can  be  determined  as 
1.8  W/cm2K.  Meanwhile,  the  maximum  heating  power  before  dry-out  occurs  is  11.6  W,  yielding  an 
interfacial  heat  flux  of  11.6/0.12=96  W/cm2.  Therefore  the  maximum  capillary  pressure  in  the  nanopores 
can  be  determined  as  -300  kPa  using  the  similar  analysis  as  on  100  pm  membranes.  The  pressure  in  the 
liquid  near  the  menisci  is  —300  kPa.  Note  that  the  overall  heat  flux  before  dry-out  is  only  -23  W/cm2  due 
to  the  porosity  of  24%.  Note  that  during  the  experiments,  the  pressure  in  the  chamber  was  kept  very  low 
(60  Pa-200  Pa)  for  efficient  removal  of  the  vapor  to  eliminate  the  thermal  resistance  associated  with 
vapor  transport.  Such  pressures  are  significantly  lower  than  the  saturation  pressure  based  on  the 
temperature  of  the  membrane.  Therefore,  the  liquid  will  evaporate  due  to  the  low  pressure  even  without 
heating,  called  “flashing”.  The  contribution  of  flashing  can  be  estimated  as  follows. 


For  100  pm  thick  membranes,  the  mass  flux  driven  by  the  pressure  difference  without  heating  will  be 

q  i-2  _ _ (21) 


m  = 


lfg 


7.33  x  105 


=  1.63  x  10-6  kg/s 


A  fraction  of  this  liquid  flux  will  be  evaporated  while  the  other  part  will  be  cooled  from  the  original 
temperature  of  the  reservoir  (22  °C)  to  the  saturation  temperature  (-20  °C)  corresponding  to  the  chamber 
pressure  (60  Pa).  The  sensible  cooling  of  the  liquid  will  provide  the  latent  heat  required  for  evaporation. 
As  a  result,  the  evaporation  mass  flux,  mevp ,  will  satisfy 

7.33  x  10 smevp  =  (22  +  20) (m  -  mevp)Cp  (22) 

where  Cp  is  the  specific  heat  of  IPA.  Therefore 
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mevp  =  2.5  x  10  ~7kg/s 

The  corresponding  heat  flux  driven  by  flashing  will  be 

q/ih  =  2.62  W/cm2 

which  is  ~  5%  of  the  total  maximum  interfacial  heat  flux  with  heating. 

For  the  50  pm  thick  membranes,  using  a  similar  calculation  method,  the  heat  flux  driven  by  flashing  will 
be  determined  as  q/th  —  5.92  W/cm2,  which  is  still  ~  5%  of  the  maximum  interfacial  heat  flux.  In 
addition  to  modeling,  a  baseline  experiment  was  performed  where  a  graduated  cylinder  containing  IPA  is 
placed  in  the  chamber  pumped  down  to  60  Pa.  The  liquid  will  evaporate  with  decreasing  rate  as  it  cools 
down  and  eventually  stops  evaporating  at  the  saturation  temperature  of  -20  °C.  At  the  beginning  of  the 
process,  the  fastest  evaporation  rate  was  found  to  be  ~  lmL/min  with  the  opening  area  of  the  cylinder 
measured  to  be  6.4  cm2.  The  evaporative  heat  flux  was  determined  as  733x0.786/60/6.4=1.5  W/cm2, 
which  is  on  the  same  order  with  model  estimations.  Therefore,  we  can  conclude  the  flashing  of  the  liquid 
is  not  significant  in  contributing  to  the  heat  transfer  of  the  total  evaporation  process. 

To  place  the  experimental  results  in  context,  we  also  developed  a  model  based  on  energy  conservation  to 
demonstrate  the  potential  of  our  proposed  configuration  to  achieve  high  heat  fluxes  shown  in  Fig.  24.  The 
permeability  of  the  membranes,  a,  can  be  determined  based  on  the  heating  power  and  driving  pressure  at 
the  transition  between  the  flooding  and  the  thin  film  regimes  in  the  experiments:  a  =  (AP /L)/(q  / phfg), 
where  A P  is  the  pressure  difference  across  the  membrane,  L  is  the  thickness  of  the  membrane,  q  is  the 
interfacial  heat  flux,  p  is  the  liquid  density,  and  hfg  is  the  latent  heat  of  the  fluid.  Therefore,  the  maximum 

interfacial  heat  flux  is  inversely  proportional  to  the  thickness  of  the  membrane,  q  =  phfgAP/aL ,  (solid 
lines  in  Fig.24).  Using  the  experimentally  determined  permeabilities  and  capillary  pressures  of  the  100 
pm  and  50  pm  thick  membranes,  an  upper  and  lower  bound  on  the  maximum  interfacial  heat  fluxes  were 
used  due  to  the  inevitable  variability  of  the  pore  geometry  during  fabrication.  The  experimental  results 
(diamond  symbols)  show  good  agreement  with  the  model.  If  membrane  thicknesses  on  the  order  of  1  pm 
were  used,  the  maximum  interfacial  heat  fluxes  could  be  as  high  as  5000  W/cm2  with  overall  heat  fluxes 
over  1000  W/cm2  even  with  membrane  porosities  of  -24%. _ 


Fig.  24  Maximum  interfacial  heat  fluxes  dependent  on  both  membrane  thickness  and  accommodation  coefficient. 
The  inverse  proportional  relation  between  heat  fluxes  and  thickness  was  validated  by  experimental  results  on  100 
and  50  pm  membranes.  By  reducing  the  thickness  and  choosing  fluid  with  higher  accommodation  coefficient, 
interfacial  heat  fluxes  as  high  as  5000  W/cm2  could  be  achieved. 
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In  previous  experiments,  the  maximum  heat  fluxes  are  only  limited  by  the  supply  of  the  fluid,  which  is 
true  in  relatively  lower  heat  fluxes.  Under  very  high  heat  fluxes,  however,  kinetic  theory  needs  to  be 
considered  in  determining  the  maximum  allowed  heat  flux  [34].  Based  on  kinetic  theory,  the  maximum 
heat  flux  across  the  interface  can  be  determined  as 


Ri  = 


2d 

2  —  d 
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fg 


M 


[pi  ,corr 


Vv 


(23) 


42tt/?  v  jfL  jtj 

where  PiiCorr  is  the  corrected  liquid  pressure  based  on  Kelvin’s  equation  where  the  curvature  of  the 
meniscus  is  accounted  for 


(  2  nvM\ 

Pi,corr -Piexpy  rpflpJ  (24) 

and  d  is  the  accommodation  coefficient  which  describes  the  probability  of  a  liquid  molecule  actually 
evaporates  in  contrast  to  reflect  back  to  liquid  phase.  The  value  of  the  accommodation  coefficient  is 
highly  sensitive  to  the  purity  of  the  system  as  well  as  the  dynamic  property  of  the  fluid.  For  polar  fluids 
such  as  ethanol,  methanol,  n-propyl  alcohol,  and  water,  the  reported  values  of  d  are  typically  0.02  to  0.04. 
For  very  pure  non-polar  fluids,  such  as  benzene  or  carbon  tetrachloride,  d  can  be  very  close  to  1  [35]. 

As  shown  with  the  membranes  tested  in  previous  experiments  (50  pm  and  100  pm  thick),  the  interfacial 
heat  fluxes  were  determined  by  the  transport  of  liquid  in  the  nanopores.  The  maximum  demonstrated 
interfacial  heat  flux  of  96  W/cm2  is  still  well  below  the  kinetic  limit  given  accommodation  coefficients  in 
the  range  from  0.02  to  0.04.  The  results  using  the  two  kinds  of  membranes  demonstrate  the  idea  that  the 
interfacial  heat  fluxes  can  be  significantly  enhanced  by  scaling  down  the  membrane  thickness. 
Considering  the  porosity  of  24%,  to  dissipate  heat  flux  over  1000  W/cm2  would  require  interfacial  heat 
fluxes  as  high  as  5000  W/cm2.  Under  such  high  heat  fluxes,  kinetic  limit  must  be  considered  where  ultra 
pure  non-polar  fluid  with  accommodation  coefficients  larger  than  0.4  will  be  required  as  well  as 
membrane  thickness  ~  1  pm.  However,  aluminium  oxide  membranes  with  thickness  smaller  than  50  pm 
are  very  flimsy  and  cannot  sustain  the  large  pressure  difference  in  the  experiments.  Future  actual  devices 
targeting  heat  fluxes  above  1000 W/cm2  would  need  supporting  structures  to  enhance  the  mechanical 
integrity  of  the  thin  membranes. 


4)  Increased  Nucleation  Density  on  Oil-Infused  Superhydrophobic  Surface  for  Enhanced 
Condensation  Heat  Transfer 

In  this  work,  we  showed  that  the  nucleation  density  on  oil-infused,  silane-coated  structured  surfaces  can 
be  significantly  increased  by  the  use  of  disordered  long-chain  silane  coatings  that  result  in  nucleation  sites 
limited  only  by  the  density  of  pillar  structures  comprising  the  surface.  We  find  that  the  increase  in  the 
nucleation  density  can  be  explained  by  heterogeneity  in  the  surface  energy  of  the  silane  coating  and  the 
reduced  water-oil  interfacial  energy.  This  effect  could  potentially  be  used  to  significantly  improve  the 
heat  transfer  coefficient  in  condensation  by  controlling  the  nucleation  density. 

Well-defined  silicon  micro/nanopillar  arrays  with  diameters,  d ,  ranging  from  0.4  pm  to  5  pm,  periods,  /, 
ranging  from  4  pm  to  25  pm,  and  heights,  h  ranging  from  10  pm  to  25  pm  were  used  in  our  experiments. 
The  silicon  surfaces  were  functionalized  with  three  different  chemicals:  1)  (Tridecafluoro-1, 1,2,2- 
tetrahydrooctyl)-l-trichlorosilane  (TFTS)  [UCT  Specialties],  which  forms  a  self-assembled  coating 
(SAC)  by  chemical  vapor  deposition  (CVD)  with  a  relatively  long  carbon  chain  (MW  =  481.54  g/mol),  2) 
Dimethyldicholorosilane  (DMCS)  [Sigma- Aldrich],  which  forms  a  self-assembled  monolayer  (SAM)  by 
CVD  with  a  short  carbon  chain  (MW  =  129.06  g/mol),  and  3)  Poly(lH,  1H,  2H,  2H-perfluorodecyl 
acrylate)  (PFDA)  polymer,  which  was  deposited  using  initiated  chemical  vapor  deposition  (iCVD)  with  a 
typical  film  thickness  of  35  nm.  Goniometric  measurements  on  smooth  functionalized  silicon  surfaces 
showed  advancing  and  receding  contact  angles  of:  9J  9V  =  122°  +  1.3778°  +  1.3°  (equilibrium  contact 
angle  9 e  ~  102.1°  ±  0.9°);  9J  9V  =  103.8°  +  0.57102.7°  +  0.4°  (0e  ~  103.2°  ±  0.3°);  and  9J  9X  =  121.1°  + 
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2.27106.3°  +  2.4°  (<9e  ~  113.5°  ±  1.6°)  for  deposited  films  of  TFTS,  DMCS,  and  PFDA  respectively.  A 
small  droplet  of  Krytox  GPL  100  oil  [DuPont]  was  applied  to  the  functionalized  silicon  pillar  arrays.  The 
surface  tension  of  Krytox  oil  is  -17-19  mN/m,  allowing  the  oil  to  spread  on  the  surface.  A  dry  nitrogen 
stream  was  used  to  assist  spreading  and  remove  excess  oil.  Typical  scanning  electron  micrographs  (SEM) 
of  the  silicon  pillar  arrays  without  and  with  the  oil  are  as  shown  in  Fig.  25(a)  and  (b). 
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Fig.  25  Scanning  Electron  Micrographs  (SEM)  of  a  pillar- stmctured  silicon  surface  (d=5  pm,  /  =  25  pm,  h  =15  pm) 
(a)  without  and  (b)  with  the  oil  infusion. _ 

We  investigated  the  nucleation  behavior  on  the  surfaces  with  and  without  the  oil  under  white  light  optical 
microscopy  (OM).  The  samples  were  horizontally  mounted  on  a  thermal  stage  [Instec  Inc.]  inside  an 
enclosure  and  cooled  to  Tw  =  283.1  +  0.1  K  in  a  dry  nitrogen  atmosphere.  Following  thermal  equilibration 
(~  5  min),  nucleation  was  initiated  by  flowing  water-saturated  nitrogen  into  the  enclosure.  The 
supersaturation,  defined  as  the  ratio  of  the  vapor  pressure  to  the  saturation  pressure  at  the  stage 
temperature  ( py/pw ),  was  controlled  by  the  temperature  of  the  water  reservoir  through  which  the  nitrogen 
carrier  gas  was  sparged  and  measured  using  a  humidity  probe  [Hygroclip,  Rotronic]  located  -1  cm  above 
the  sample.  Typical  values  of  the  supersaturation  were  around  S  ~  1.6.  The  nucleation  density  and 
subsequent  growth  behavior  was  recorded  at  a  frame  rate  of  30  fps  using  a  CMOS  camera  [Phantom 
V7.1,  Vision  Research]  attached  to  the  optical  microscope.  As  shown  in  Fig.  26  (a)  and  (b),  the  nucleation 
density  on  oil-infused  surface  with  long-chain  TFTS  coating  was  significantly  higher  than  that  on  the 
same  surface  without  oil.  We  did  not  observe  any  nucleation  between  pillars  which  can  be  attributed  to 
the  large  thickness  of  oil  coverage  between  pillars  that  introduces  a  large  diffusion  resistance  for  water 
vapor.  However,  a  droplet  was  formed  on  the  tip  of  almost  each  pillar  where  the  oil  film  was  thin  enough 
to  allow  sufficient  mass  to  diffuse  to  the  functionalized  silicon  surface.  Meanwhile,  on  pillar  arrays  coated 
with  DMCS  and  PFDA,  we  did  not  observe  similar  increases  in  nucleation  density,  as  shown  in 
Fig.  26(c). 
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Figure  26(a)  White-light  microscopy  image  of  condensation  on  a  silicon  pillar  array  without  oil,  with  nucleation 
sites  highlighted  by  red  circles,  (b)  White-light  microscopy  image  of  condensation  on  a  composite  surface  with  oil. 
Nucleation  occurred  on  the  tip  of  almost  every  pillar,  (c)  Comparison  of  the  change  in  nucleation  density  with 
various  surface  coatings.  The  nucleation  densities  were  normalized  against  pillar  densities  for  a  fair  comparison. 


To  investigate  the  mechanism  for  this  drastic  change  in  nucleation  density,  we  performed  atomic  force 
microscopy  (AFM)  in  tapping  mode  on  a  smooth  TFTS -coated  silicon  surface  and  observed  the  presence 
of  nanoscale  agglomeration  of  silane  molecules,  as  shown  in  Fig.  27(a).  Such  agglomerates  have  been 
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observed  in  previous  studies  and  considered  as  disordered  agglomeration  of  excessive  silane  molecules 
[36].  The  phase  image  of  the  AFM  measurement  (Fig.  27(d))  showed  significantly  higher  phase  angle  on 
the  agglomerates,  which  suggests  that  the  agglomerates  were  locally  more  hydrophilic  compared  to  the 
background  film.  Such  heterogeneity  was  also  supported  by  the  high  contact  angle  hysteresis  (0a  /  0T  = 
122°  +  1.3778°  +  1.3°).  However,  such  agglomerates  were  not  observed  on  other  kinds  of  surface 
coatings  such  as  DMCS  and  PFDA,  as  shown  in  Fig.  27(b)  and  (c). 
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Fig.  27  Atomic  Force  Microscopy  (AFM)  height  images  of  smooth  silicon  surfaces  coated  with  (a)  TFTS,  (b) 
DMCS,  and  (c)  PFDA.  Micelles  structures  were  only  observed  on  TFTS  coated  surfaces,  (b)  AFM  phase  image  of 
smooth  silicon  surface  coated  with  TFTS.  The  local  high  phase  angle  indicates  higher  hydrophilicity  of 
agglomerates. 


The  local  contact  angles  on  the  hydrophobic  substrate  and  the  nanoscale  agglomerates  can  be  determined 
based  on  a  modified  Cassie-Baxter  model.  Assuming  the  local  contact  angles  on  the  hydrophobic 
substrate  and  the  nanoscale  agglomerates  to  be  0?  and  82,  respectively,  the  macroscopic  advancing  and 
receding  contact  angles  are  determined  as 

9a  =  0i  (25) 

cos8r  =  Jfcos62  +  (l  —  sff  )cos01  (26) 

where/ is  the  area  fraction  of  the  agglomerates  [37]. 

Based  on  the  macroscopically  measured  advancing  and  receding  angles,  0a  =  122°  +  1.3°  and  6X  =  78°  + 
1.3°,  and  the  fraction  of  the  agglomerates  determined  as /~  0.4  from  AFM,  the  local  contact  angles  on  the 
hydrophobic  substrate  and  the  hydrophilic  micelles  were  found  to  be  9\  =  122°  +  1.3°  and  02  =  60°  +  1.5°. 

The  nucleation  rate,  7,  can  be  determined  by  classical  nucleation  theory  (CNT)  as  [38] 
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(27) 


/  =  zf*  exp(-G*) 

In  Eqn.  (24),  z  is  the  Zeldovich  factor  and  G*  is  the  dimensionless  energy  barrier,  given  by 

z  =  ( kT  In  S)2/8nv0J  kTip(6)y3 

G*  =  167nK0KV /3(/c7)2  (lnS)2 

where  S  is  the  super  saturation  and  is  the  activity  that  accounts  for  the  effect  of  contact  angle,  f*  is 
the  frequency  of  monomer  attachment  to  the  critical  droplet  nucleus  dependent  on  the  nature  of  the 
nucleus  growth.  The  main  modes  of  growth  during  heterogeneous  nucleation  are  limited  via  surface 
diffusion  or  direct  impingement  of  monomers  to  the  nucleus  [39,  40].  Volumetric  diffusion  is  a  third 
growth  limiting  step,  which  is  only  considered  important  for  nucleation  taking  place  in  liquid  or  solid 
solutions  [6].  However,  we  included  all  three  mechanisms  when  calculating  the  nucleation  rates. 

The  frequency  of  monomer  attachment  due  to  direct  vapor  impingement  is  given  by 

f\  =  7n[(l  “  cos(0w))/2t/>2/3(0)] (36uVq  )1/3/n2/3  (30) 

where  yn  is  the  sticking  coefficient  (0  <  yn  <  1),  I  is  the  classical  Hertz-Knudsen  impingement  rate 
(/  =  P / \[2nm0kT  )),  n  is  the  number  of  molecules  in  the  nucleated  cluster,  and  v0  is  the  volume  of  an 
individual  water  molecule  (v0  =  3  x  10"29  m3).  To  determine  an  upper  bound  on  the  nucleation  rate,  a 
sticking  coefficient  of  one  was  assumed  (yn  =  1). 

The  frequency  of  monomer  attachment  due  to  surface  diffusion  is  given  by  /*sd  =  ync*As  /,  where  c*  is 
the  capture  number  due  to  surface  diffusion  (1  <  c*  <  5),  and  2S  is  the  mean  surface  diffusion  distance  of 
an  adsorbed  monomer  on  the  substrate.  The  capture  number  c*  is  size  independent  and  approximately 
equal  to  1.9  for  heterogeneous  condensation  of  water  vapor  [41].  The  mean  surface  diffusion  distance  is 
dependent  on  the  wettability  of  the  substrate  and  is  given  by  As  =  where  Ad  is  the  surface 

diffusion  coefficient  (Dsd  =  ds  vs  exp[—Esd/kT]),  rd  is  the  desorption  time  (rd  =  (1/  vs)  exp[— Edes/ 
kT ]),  vs  is  the  adsorbed  molecule  vibration  frequency  determined  using  the  Debye  approximation 
(vs  =  VDa /  2),  ds  is  the  length  of  a  molecular  jump  along  the  substrate  surface  approximated  by  the 
lattice  constant  of  the  substrate  ( ds  =  5.4  A)  [42]  and  VD  is  the  speed  of  sound  in  the  substrate  (Vb  =  8433 
m/s).  The  desorption  and  surface  diffusion  energies  are  given  by  Edes  =  E1  +  <rsva0  and  Ad  —  0.5Aies> 
respectively,  where  Ej  is  the  binding  energy  of  an  n  =  1  sized  cluster,  crsv  is  the  solid  vapour  interfacial 
energy  and  a0  is  the  water  molecule  surface  area  ( a0  =  4.67  x  10"19  m2).  The  calculated  energies  of 
desorption  show  excellent  agreement  with  that  of  experiment  and  molecular  dynamics  simulations 
(Aes,si02=  0.9  eV)  [43]. 


(28) 

(29) 


The  frequency  of  monomer  attachment  due  to  volumetric  diffusion  is  given  by 

f\ d  =  Yn  (l  -  C;S1703W)  (6?r2v oV^DCn1/3  (31) 

where  D  is  the  self  diffusion  coefficient  of  water  vapor  ( D  =  (3/  QnriQd^^kT / nm0 ),  C  is  the 
equilibrium  concentration  of  monomers  ( C  =  (l/a0)  exp(— W1/kT)),  d0 ,  m0  and  n0  are  the  water 
molecule  diameter  (d0  =  3.0  A)  [43],  mass  (m0  =  3  x  10"26  kg)  and  number  density  (n0  =  Aa/vm), 
respectively. 

By  adding  the  nucleation  rate  from  the  three  mechanisms  together,  the  nucleation  rate,  /,  can  be 
determined  as  a  function  of  the  contact  angle  and  surface  tension  of  the  condensate  at  given 


Approved  for  public  release:  distribution  unlimited 


28 


supersaturations,  as  shown  in  Fig.  28. 


Fig.  28  Nucleation  rate  predicted  by  the  classical  nucleation  theory  as  a  function  of  contact  angle  and  interfacial 
tension  at  S  =  1.6. 


The  surface  tension  of  water  in  air  is  72  mN/m  and  the  interfacial  tension  between  water  and  the  Krytox 
oil  was  found  to  be  58  mN/m  by  measuring  the  contact  angle  of  water  droplet  on  oil  film.  From  Fig.  28 
we  can  see  that  for  the  hydrophilic  micelles  ( 9  ~  60°),  the  reduced  interfacial  tension  between  water  and 
oil  leads  to  a  significant  increase  in  the  nucleation  rate.  With  the  nanoscale  agglomerates  acting  as 
nucleation  sites,  nucleation  was  supposed  to  occur  on  almost  every  pillar  tip  where  the  oil  film  was  thin 
enough  for  water  vapor  to  diffuse  through.  On  surfaces  without  the  agglomerates s,  such  as  DMCS  and 
PFDA-coated  surfaces,  the  contact  angles  are  over  100°  and  the  nucleation  rate  was  essentially  zero  even 
with  reduced  interfacial  tension. 


(a)  (b) 


Fig.  29  Optical  microscope  images  of  the  condensation  experiments  on  3-TMPM  coated  silicon  pillar  arrays  (a) 
without  (N  =  6.7xl08m'2)  and  (b)  with  Krytox  oil  (N  =  1.7xl010m'2).  Nuclei  are  highlighted  by  red  circles  in  (a)  and 
pillars  without  nucleation  were  highlighted  by  blue  circles  in  (b).  The  pillars  have  diameters  of  2.5  pm,  periods  of 
7.5  pm  and  heights  of  25  pm.  (T2  =1.7xlQ10m~2). 


In  order  to  validate  our  assumption,  we  carried  out  condensation  experiments  using  silicon  pillar  arrays 
coated  with  3-(trimethoxysilyl)propyl  methacrylate  (3-TMPM).  The  advancing  and  receding  contact 
angle  of  water  on  a  smooth  silicon  surface  coated  with  3-TMPM  are  65°  +  1.5°  and  53°+  1.1°, 
respectively.  The  contact  angle  is  in  the  range  where  the  nucleation  rate  will  be  almost  zero  with  a  surface 
tension  of  72  mN/m  and  nucleation  should  occur  on  every  tip  of  pillars  with  an  interfacial  tension  of  58 
mN/m.  Optical  images  of  the  condensation  experiments  are  as  shown  in  Fig.  29.  Similar  to  the  behavior 
of  TFTS-coated  surface,  a  significant  increase  in  nucleation  density  was  observed  with  the  addition  of 
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Krytox  oil  as  expected.  This  result  supports  well  our  hypothesis  related  to  the  role  of  the  hydrophilic 
nanoscale  agglomerates  in  the  droplet  nucleation  process. 

Next,  we  studied  the  overall  heat  transfer  performance  of  an  immersion  condensation  surface.  While 
studies  on  well-defined  silicon  micropillar  arrays  can  provide  physical  insight  into  immersion 
condensation  behavior,  they  are  not  practical  due  to  cost  and  challenges  in  interfacing  the  silicon  substrate 
and  our  heat  transfer  measurement  apparatus  with  minimum  uncertainties.  Therefore,  we  performed 
immersion  condensation  heat  transfer  measurements  on  oil-infused  copper  oxide  (CuO)  nanostructures 
functionalized  with  TFTS,  which  promises  a  scalable,  low  cost  platform  for  condensation  surfaces.  [30] 
SEM  images  of  representative  copper  oxide  nanostructures  without  and  with  Krytox  oil-infusion  are 
shown  in  Fig.  30a  &  b,  respectively.  Condensation  experiments  were  performed  on  the  CuO  surfaces 
without  and  with  oil-infusion  in  an  environmental  SEM  with  1  <  S  <  1.29  for  visualization.  [12]  The 
Figures  30c  &  d  show  an  order  of  magnitude  increase  in  nucleation  density  on  the  oil-infused  surface,  as 
similarly  observed  on  the  silicon-based  microstructures.  To  capture  the  condensation  heat  transfer 
behavior,  we  formed  the  oil-infused  heterogeneous  CuO  surfaces  on  copper  tubes  (see  Methods  for 
detailed  fabrication  process).  Figures  30e  &  f  show  condensation  on  a  typical  hydrophobic  surface  and  an 
oil-infused  heterogeneous  surface,  respectively.  Significantly  higher  droplet  densities  were  observed  on 
the  oil-infused  surface,  while  a  low  departure  radius  of  approximately  0.5  mm  was  maintained.  Prior  to 
droplet  departure,  the  droplets  grew  orders  of  magnitude  larger  than  the  characteristic  length  scale  of  the 
nanostructures,  thus  high  apparent  contact  angles  of  the  droplet  (-110°)  were  observed,  consistent  with 
the  low  surface  energy  of  the  solid-oil  composite  surface. 


Fig.  30  Scalable  copper  oxide  (CuO)  surfaces  for  immersion  condensation,  (a)  Field  emission  SEM  (FESEM)  image 
of  CuO  nanostructures,  (b)  Environmental  SEM  (ESEM)  image  of  CuO  nanostructures  infused  with  Krytox  oil.  (c) 
ESEM  image  of  nucleation  on  TFTS-coated  CuO  surface,  (d)  ESEM  image  of  oil-infused  TFTS-coated  CuO 
surface.  An  order  of  magnitude  higher  nucleation  density  was  observed  compared  to  (c).  (e)  Image  of  drop  wise 
condensation  on  a  hydrophobic  copper  tube  surface,  (f)  Image  of  condensation  on  an  oil -infused  TFTS-coated  CuO 
surface.  Significantly  higher  droplet  density  was  observed  on  the  oil-infused  surface  while  a  low  departure  radius  of 
approximately  0.5  mm  was  maintained. 


Overall  heat  transfer  coefficients  were  measured  to  evaluate  the  performance  on  three  different  CuO- 
based  surfaces:  a  hydrophobic  surface  for  typical  drop  wise  condensation  (DHP),  a  superhydrophobic 


Approved  for  public  release:  distribution  unlimited 


30 


TFTS-coated  copper  oxide  surface  (SHP),  and  a  Krytox  oil-infused,  TFTS-coated  CuO  surface  (Fig.  31). 
The  Krytox  GPL  100  oil  evaporates  completely  when  the  test  chamber  is  evacuated  to  pressures  lower 
than  1  Pa.  Therefore,  we  set  the  initial  chamber  pressure  as  high  as  30  Pa  (primarily  composed  of  non- 
condensable  gases,  NCG)  to  avoid  the  evaporation  of  oil  with  steam  pressures  ranging  from  2  to  3  kPa  (1 
<  S  <  1.6)  in  the  experiments.  This  is  consistent  with  actual  condenser  systems  where  NCG  partial 
pressures  are  typically  found  in  the  range  of  30  Pa  and  significantly  affect  the  condensation  heat  transfer 
performance.  [44-47]  Accordingly,  with  these  experimental  conditions,  we  were  able  to  emulate  a  more 
realistic  condensation  environment  and  demonstrate  the  practical  significance  of  the  immersion 
condensation  mode.  While  the  SHP  surface  is  more  hydrophobic  than  the  DHP  surface,  flooding  and 
strong  pinning  of  the  condensate  was  observed  due  to  the  high  supersaturation  conditions  ( S  as  high  as 
1.6),  leading  to  similar  heat  transfer  coefficients  with  the  DHP  surfaces.  Note  that  these  results  are 
distinct  from  previous  literature  where  jumping  of  droplets  on  SHP  surfaces  increased  heat  transfer 
coefficients  at  lower  saturation  conditions  ( S  <  1.1 2).  [48]  In  addition,  the  overall  heat  transfer  coefficients 
on  DHP  surfaces  in  this  work  (/z~2-7  kW/m2K)  are  much  lower  compared  to  pure  vapor  conditions  (h~  12- 
13  kW/m2K)[48]  due  to  the  presence  of  NCGs  acting  as  a  diffusion  barrier  to  the  transport  of  water  vapor 
towards  the  condensing  surface.  In  comparison  to  the  DHP  surfaces,  the  Krytox  oil-infused  TFTS-coated 
CuO  surface  demonstrated  approximately  a  100%  improvement  in  heat  transfer  coefficient  over  the  entire 
range  of  supersaturations  tested  (1  <  S  <  1.6)  with  the  existence  of  NCGs.  While  the  available 
condensation  area  was  reduced  due  to  the  significant  oil  coverage,  the  significant  improvement  in  the 
overall  heat  transfer  coefficient  highlights  the  collective  role  of  enhanced  nucleation  density,  more 
frequent  droplet  removal,  and  lower  droplet  contact  angle. _ 


Fig.  3 1  Experimental  immersion  condensation  heat  transfer  measurement.  Comparison  of  overall  heat  transfer 
coefficient  during  condensation  on  the  hydrophobic  surface,  TFTS-coated  superhydrophobic  surface,  and  oil-infused 
composite  surface  with  30  Pa  of  non-condensable  gases.  The  supersaturation  was  varied  in  the  range  1<S<1.6.  The 
heat  transfer  coefficient  on  the  oil-infused  surface  increased  by  approximately  100%  compared  to  the  drop  wise  and 
superhydrophobic  surfaces. 


Oil-impregnated  surfaces  have  been  recently  reported  as  a  promising  approach  to  enhance  condensation 
heat  transfer  surfaces  due  to  the  ultra-low  droplet  adhesion.  [49-51]  However,  easy  droplet  removal  is  not 
the  only  desired  property  for  high  heat  transfer  performance.  Low  contact  angle  and  high  nucleation 
densities  are  also  essential  for  maximized  condensation  heat  transfer.  In  this  work,  we  have  demonstrated 
that  by  combining  surface  heterogeneity  and  oil-infusion,  the  nucleation  density  in  condensation  can  be 
increased  by  over  an  order  of  magnitude  via  immersion  condensation  while  maintaining  low  droplet 
adhesion,  which  has  never  been  observed  previously.  The  increase  in  nucleation  densities  was  achieved 
by  a  combined  effect  of  heterogeneity  and  reduced  oil-water  interfacial  tension  as  explained  by  our  model 
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based  on  classical  nucleation  theory  and  the  model  was  validated  by  control  experiments  using  silane- 
coated  silicon  micropillar  arrays.  With  improved  understanding  of  the  physics,  we  investigated  oil-infused 
superhydrophobic  copper  oxide  surfaces  as  a  platform  for  condensation  enhancement  in  practical  systems. 
We  demonstrated  that  the  condensation  heat  transfer  coefficient  on  such  oil-infused  heterogeneous 
surfaces  can  be  enhanced  by  approximately  100%  compared  to  state-of-the-art  drop  wise  surfaces  in  the 
presence  of  non-condensables  gases.  Further  work  is  needed  to  tailor  oil  and  coating  properties,  as  well  as 
surface  geometry  to  minimize  oil  loss  during  operation  and  maximize  condensing  surface  area.  With 
continued  development,  immersion  condensation  promises  to  be  an  important  condensation  mode  for  a 
variety  of  heat  transfer  and  resource  conserving  applications. 
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